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Abstract 



We consider a coupled system of two singularly perturbed reaction-diffusion equations, with 
two small parameters < e < fi < 1, each multiplying the highest derivative in the equations. 
The presence of these parameters causes the solution(s) to have boundary layers which overlap 
and interact, based on the relative size of e and fi. We construct full asymptotic expansions 
together with error bounds that cover the complete range < e < [i < 1. For the present case 
of analytic input data, we derive derivative growth estimates for the terms of the asymptotic 
expansion that are explicit in the perturbation parameters and the expansion order. 



1 Introduction 



Singularly perturbed (SP) boundary value problems (BVPs), and their numerical approximation, 
have received a lot of attention in the last few decades (see, e.g., the classical texts |13|. [T7j on 
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asymptotic analysis and the books [TO], [TT], [TJ], whose focus is more on numerical methods for 
this problem class). One common feature that these problems share is the presence of boundary 
layers in the solution. In order for a numerical method, designed for the approximation of the 
solution to SP BVPs, to be considered robust it must be able to perform well independently of the 
singular perturbation parameter (s). To achieve this, information about the regularity of the exact 
solution is utilized, and in particular, bounds on the derivatives. Such information is available 
in the literature for scalar SP BPVs of reaction- and convection-diffusion type in one- and two- 
dimensions (see, e.g., [SI [7] for scalar versions of the problem studied in the present article). For 
systems of SP BVPs, the bibliography is scarce, even in one-dimension; see the relatively recent 
review article [4] and the references therein, for such results available to date. It is, therefore, the 
purpose of this article to begin filling this void; in particular, we provide the regularity theory 
for a system of two coupled SP linear reaction-diffusion equations, with two singular perturbation 
parameters. Our analysis is complete for the problem under consideration in that we derive full 
asymptotic expansions for all relevant cases of singular perturbation parameters and give explicit 
control of all derivatives of all terms appearing in the expansions. Even though this is a linear, 
one-dimensional problem, the methodology presented here can be the starting point for treating 
more difficult problems. 

The regularity results obtained here are used in the companion communication [9j to prove, for the 
first time, exponential convergence of the /ip-FEM for problems with multiple singular perturbation 
parameters. This exponential convergence result for the /ip-FEM relies on mesh design principles 
firmly established for problems with a single singular perturbation parameter as discussed in \15\ 
EEl El [SI [6] ; the mathematical analysis of [9] shows that these mesh design principles extend to 
problems with multiple singular perturbation parameters and confirms the numerical results of [18] 
for the problem class under consideration here. 

The rest of the paper is organized as follows: In Section[2]we present the model problem and discuss 
the typical phenomena. In Sections [SHE] we address the regularity of the solution, as it depends on 
the relationship between the singular perturbation parameters. The proofs of most of the results 
presented in these sections are technical and are relegated to Appendices lATTDl 

In what follows, the space of square integrable functions on an interval I C R will be denoted by 
I? (I), with associated inner product 



We will also utilize the usual Sobolev space notation H k (I) to denote the space of functions on 



where dl denotes the boundary of /. For z £ C, we use dB r (z) to denote the ball of radius r 
centered at z. Finally, the letter C will be used to denote a generic positive constant, independent 
of any discretization or singular perturbation parameters and possibly having different values in 
each occurrence. 





itf (!) = {«€ IT 1 (I):u\ dI = 0}, 
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2 The Model Problem and Main Results 



We consider the following model problem: Find a pair of functions (u, v) T such that 

-e 2 u"(x) + axi(x)u(x) + a\2{x)v(x) = f{x) in / = (0, 1), 
-fi 2 v"(x) + a%\{x)u{x) + a,22(x)v(x) = g(x) in I = (0, 1), 



along with the boundary conditions 

w(0) = u(l) = , v{0) = v(l) 

With the abbreviations 

s 2 



0. 



(2.1a) 



(2.1b) 



U 



ll 2 



A(x) 



an (a;) 012(2;) 
a 21 (x) a 2 2(x) 



equations (I2.1ap - (|2.1bj) may also be written in the following, more compact form: 

-E £ '^XJ"(x) + A(x)U = F, U(0) = U(l) = 0. 



(2.2) 



The parameters < e < \i < 1 are given, as are the functions /, g, and Oy, i,j G {1,2}, 
which are assumed to be analytic on I = [0, 1]. Moreover we assume that there exist constants 
C f> 7/, C g , 7 9 , C a , 7 a > such that 



» 



/ (n) || Loo(/) <C fl n f n\ V n G No, 
9 {n) \\ L ^ {I) <C gl n g n\ VneNo, 

<C a7 >! VnGN ,*,jG{l,2} 



(2.3) 



L°°(J) 



The variational formulation of (|2.1ap - (|2.1bp reads: Find U := 



[u,v 



G [tfo 1 (I)] 2 such that 



(2.4) 



B(U,V)=F(V) V V : = («, t?) G (/)] 2 , 

where, with (•, the usual L 2 (I) inner product, 

2?(U,V) = e 2 (u' ,u) j + fi 2 ,v') j + (anu + ai 2 v,u)j + (a2iu + a 2 2V,v)j , (2.5) 
F(V) = (f^j + ig^j. (2.6) 

The matrix- valued function A is assumed to be pointwise positive definite, i.e., for some fixed a > 

~t T A~t > a 2 ~i T ~i V^eR 2 Vx G 7. (2.7) 

It follows that the bilinear form B (•, •) given by (|2.5p is coercive with respect to the energy norm 



IUI 



in, u) 



l EJ 



e \u\ 



l«llo,/ + IMIo,// • 



(2.8) 



i.e., 



5(U,U)>||U||| i7 VUG [Ho 1 (J)] . 

This, along with the continuity of B (■,■) and i 7 (•) , imply the unique solvability of (|2.4p . We also 
have, by the Lax-Milgram lemma, the following a priori estimate 



|U||_g j < max <, 1, 



a 



2 + llnll 2 
0,7" + llffllo,!- 



(2.9) 
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For the development of certain asymptotic expansions, it will be convenient to observe that our 
assumption (|2.7p implies (see Lemma IA.2I for the proof ) 

a k k(x) > aVxGl, k = 1,2, (2.10) 
det A(x) = an(x)ci22(x) — ai2(x)a,2i(x) > a max{an(i), 022(2;)} > a 2 Vx G /. (2-11) 

We note that the special case when the parameters are equal, i.e. e = \i, was analyzed in [TU]. In 
the general case considered here, there are three scales (1 > fj, > e) and the regularity depends on 
how the scales are separated. Correspondingly, there are 4 cases: 



(I) The "no scale separation case" which occurs when neither \xf \ nor e/fi is small. 

(II) The "3-scale case" in which all scales are separated and occurs when \xj\ is small and ej ^ is 
small. 

(III) The first "2-scale case" which occurs when u/1 is not small and ej [i is small. 

(IV) The second "2-scale case" which occurs when u/1 is small and ej \i is not small. 



The concept of "small " (or "not small") mentioned above, is tied in two ways to our performing 
regularity theory in terms of asymptotic expansions. First, on the level of constructing asymptotic 
expansions, the decision which parameters are deemed small determines the ansatz to be made 
and thus the form of the expansion. Second, on the level of applying asymptotic expansions, the 
decision which parameters are deemed small depends on whether the remainder resulting from the 
asymptotic expansion can be regarded as small. 

We need to introduce some notation: 

Definition 2.1. 1. We say that a function w is analytic with length scale v > (and analyticity 
parameters C w , ^ w ), abbreviated w G A(fi, C w , 7^), if 

||^ (n) < C wl l max{n, i/" 1 }" Vn G N . 

2. We say that that an entire function w is of L°° -boundary layer type with length scale v > 
(and analyticity parameters C w , ^ w ), abbreviated w G BC°°{y, C w , j w ), if for all x £ I 



■w 



(n) 



(x)\ < C w ~/y- n e- dist ( x ' 9I V u Vn G N . 



3. We say that that an entire function w is of L 2 -boundary layer type with length scale v > 
(and analyticity parameters C w , j w ), abbreviated w G B£ 2 (is,C w ,j w ), if for all x £ I 

\\ e ^t(x,dl)/u w (n) < (j^l^y-n y n £ Nq> 

All three definitions extend naturally to vector-valued functions by requiring the above bounds com- 
ponentwise. 

In view of the length of the article, we collect the main result at this point; the four scale separation 
cases (jl|)- pvp listed above correspond to the four cases listed in the following theorem: 
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Theorem 2.2. There exist constants C, b, 8, q, 7 > independent o/0<e< / u<l such that the 
following assertions are true: 

(I) Ue^Cr 1 /^). 

(II) U can be written as U = W+U BL +U BL +R, where W G .4(1,(7,7), U BL G B£°°(6fj,, C, 7) 
U BL G S£°°(fe,(7,7), and ||R|| LO o (97) + ||R|b,/ < C [e~ b ^ + e "^/ £ ] . Additionally, the 
second component vofXJ satisfies the sharper estimate v G BC°°(8e, C(e/n) 2 ,j). 

(Ill) If e/ ijl < q then U can be written as U = W + XJbl + R; where W G 4(^, C, 7), G 
B£°°(8e, C, 7), and [|R[| i0 o(gn + HRHe,/ < Ce~ b / £ . Additionally, the second component v of 

U satisfies the sharper estimate v G BC°°(8e,C(e/ fi) 2 ,7). 

fIVj U can be written as U = W+U BL +R, where W G 4(1, C, 7), U BL G BC°°{8pL, C^/Jxje, 7/x/e), 
and ||R|| £ oo( 8 j) + ||R||b,/ < C(n/ef e~ b ^ . 

Proof. This result is obtained by combining Theorems 13, H 14. H 15, H and 16,11 to be found in Sec- 
tions [3H3 We emphasize that some results in these theorems are slightly sharper since they analyze 
all terms of the asymptotic expansions whereas Theorem 12,21 is obtained from the asymptotic ex- 
pansions by suitably selecting the expansion order. □ 



3 The no scale separation case: Case [I] 



In this case neither fi/1 nor e/fi is small, which means that the boundary layer effects are not very 
pronounced. By the analyticity of aij,f and g, we have that u and v are analytic. Moreover, we 
have the following theorem. 

Theorem 3.1. Let U = (u,v) T be the solution to (|2.1a|) - (|2.1b|) with < e < fj, < 1. Then there 
exist constants C and K > 0, independent of e and fi, such that U G A(e, Ce~ l l 2 , K) satisfying 
the sharper estimate 

lit") + v {n) <Crmaxkr 1 }' i VnGN . (3.1) 
0,/ 0,/ 

Proof. The L 2 -based estimate (|3.ip was shown in |19j for the special case e = fi. The extension to 
current situation e < /i is straight forward. We note that the Sobolev embedding theorem in the 
form — C|MIl 2 (7) IMIh 1 ^) allows us to infer from (|3.ip the assertion U G A(e, Ce~ x l 2 , 7) 

for suitable C, 7 > independent of e and /x. □ 



4 The three scale case: Case II 



In this case all scales are separated and it occurs when both /i/1 and ej [i are deemed small. This is 
arguably the most interesting (and challenging from the approximation point of view) case, since 
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boundary layers of multiple scales appear. Additionally, this case shows most clearly the general 
procedure for obtaining asymptotic expansions and error bounds for problems with multiple scales. 
Before developing the asymptotic expansion, we formulate the main result: 



Theorem 4.1. The solution U of \2. 1\) can be written as U = W + XJbl + Ubl + R ; where 
W <E A(l,Cw,Jw), U G BC°°(5fi, Cbl, Ibl), U G B£.°°(5e, Cbl,1bl) for suitable constants Cw, 
Cbl, 7w '7 1BL> S > independent of fi and e. Furthermore, R satisfies 

||R|U~ ( aj) + ||i^R|U^ < c[e- b ^ + e- b ^ £ 
for some constants C , b > independent of /x and e. In particular, \\R\\e,i < C [e~ b ^ + e _b ^/ £ ] . 
Additionally, the second component v of TJ bl satisfies the sharper regularity assertion 

veBC°°{Se,CBL{£ln?,lBLy 

Proof. See Section 14.61 □ 

Anticipating that boundary layers of length scales O(fi) and 0(e) will appear at the endpoints 
x = and x = 1, we introduce the stretched variables x = x/fj, ,x = x/e for the expected layers at 
the left endpoint x = and variables 2^ ' = (1 — x)/(jl , x R = (1 — x)/e for the expected behavior 
at right endpoint x = 1. We make the following formal ansatz for the solution U: 

00 00 ' / \ j 

U ~ E E (I)' ( £ ) [ U ^(^) + + U^-(x) + V*(x R ) + \J*(x R )] , (4.1) 

i=0 i=0 

where the functions Uy, U^-, U^, will be determined shortly. The decomposition of 
Theorem 14.11 is obtained by truncating the asymptotic expansion (14. lh after a finite number of 
terms: 

tM/„\ ._ ( u ( x ) \ _ w („\ 1 ttM f~\ , vM /-_R\ , frAf f~\ , i>M /-~.R\ 



(4.3) 



where 

Mi M 2 

wm(*) = z;z> 

i=0 j=0 

denotes the outer (smooth) expansion, 

-££>(£)'( II ) ■ n^) ) • («> 

i=0 3=0 V ^ y v ,JV ' 7 i=0 j=0 v ^ y v l -? v y 7 

denote the left and right inner (boundary layer) expansions associated with the variables x, x R , 
respectively, 





Uij 


(x) \ 




V % 


(x) ; 




Mi 


M 2 
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denote the left and right inner (boundary layer) expansions associated with the variables x, x R 
respectively, and 



R M (x) :-- 



r u {x) 
r v (x) 



U(x) - (W M (x) +Uf L (x) + V£ L (S") + U£ L (2f) + V^ L (x«)) (4.6) 



M i^R\ 



tM 



tM i~R- 



denotes the remainder. Theorem 14.11 will be established by selecting M\ = 0(1/ ji), M2 = 0(/x/e). 



4.1 Derivation of the asymptotic expansion 



In order to derive equations for the functions U fj , U5, the procedure is as follows: 

first, the ansatz (|4.ip is inserted in the differential equation (|2.1a|) . then the scales are separated 
and finally recursions are obtained by equating like powers of fi and e/fj,. 

In order to perform the scale separation, we need to write the differential operator L e ^ in different 
ways on the various scales. In particular, for the x and the x-scales, the coefficient A is written, 
by Taylor expansion, as 



A(x) = 5> fc A fc 5* A fe :=A( fe )(0) 

k=0 

00 / v k 

fro W 



«ff(0) < 2 >(0) 
k\ k\ 

k\ fc! 



(4.7) 
(4.8) 



Corresponding representations are obtained for the variables x R and x R by expanding around the 
right endpoint x = 1. Hence, the differential operator L £j „ applied to a function depending on x or 
x takes the following form: 



on the x-scale: 



on the x-scale: 



-fi- 2 B^dl\J(x) + Y, H k A k x k XJ(x) , (4.9) 

k=0 

00 / \ k 

-e- 2 E £ '^|U(x) + ^/i fc (-) A fe x fc U(x). (4.10) 
k=o ^ 



Clearly, analogous forms exist for the operator on the x R and x R scales. We now insert the ansatz 
(|4.ip in the differential equation equation (|2.1a|) . where the differential operator L e ^ takes the form 
given above on the fast scales X , X • X ■ X , and we separate the scales, i.e., we view the variables 
X X • X ■ X X 3jS independent variables. Then, we obtain 



00 00 / \ 3 

EE"* ~ [-E^+A(x)U, 
i= oi=o v/x/ 



EE 

co 00 

EE 

j=0 j=0 



-/x- 2 E^(Ug)" + E^ A ^ fcU 

k=0 

00 



k=0 



0, 

0. 



(4.11) 
(4.12) 
(4.13) 
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and two additional equations for U^, corresponding to the scales x R , x R that are completely 
analogous to (|4TT2l) . (|4TT3l) . We write 

M:;)' <*-(?)• ( " 4) 

and equate like powers of and e/// in (14.1ip . (14.12p . fl4. 13f) to get the following recursions: 

^- 2 )+A(x)U ij = F{j, (4.15a) 

V i~2,j J 

iU k~\i" )+E A ^ fc Uf_ fcj = 0, (4.15b) 
^ Ui d> / fc=0 

^ V + E A fc x fc Uf_ fcj _ fc = 0, (4.15c) 

where we adopt the convention that if a function appears with a negative subscript, then it is 
assumed to be zero. Furthermore, we set 



' oo 



f g ) . F„ - -0 if ( i . j ) : (0.0). 



The procedure so far leads to a formal solution U of the differential equation (|2.1ap ; further 
boundary conditions are imposed in order to conform to the boundary conditions (|2.1b|) . namely, 

Uij-(O) + U§(0) + U§(0) = 0, plus decay conditions for Ujjj, Ug at +oo (4.15d) 
with analogous conditions at the right endpoint x = 1, which couple Vij, V R , and V R . 



4.2 Analysis of the functions Uj 



Since the matrix A(x) is invertible for every x £ I, equation (|4.15ap may be solved for any i, j 
yielding 



^0,0 


= A- 1 


" / " 


. ^0,0 _ 




. 9 _ 



and for / (0,0) 



with, as mentioned above, 



Uj 



*i-2j-2 
u" 

V i-%3 



(4.16) 
(4.17) 



= , = if i < or j < 0. 

Note that (|4.17p gives all the cases (£, 0) and (0,j) because the right-hand side in (|4.17p is known. 
Moreover, for each j, (I4.17P allows us to compute Uij, Vij V i, thus (|4.16l) - (l4.17p uniquely determine 

Uij , Vij V i , j. 



We have the following lemma concerning the regularity of the functions Uj,-: 
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Lemma 4.2. Let f , g, and A satisfy \2. 3\) and \2.1^ . Let Uij, Vij be the solutions of (|4.16j) . (|4.17p . 
Then there exist positive constants Cs and K and a complex neighborhood G of the closed interval 
I independent of i and j such that 



OVj >i, 



V%j — Vij 



if i or j is odd, 



Wij{z)\ + \vij{z)\ < CsS^KH* V zeG s :={zeG: dist(z,dG) > 5}. 



(4.18) 

(4.19) 
(4.20) 



Proof. The proof is by induction on i, where the estimate fj4.20f) follows by arguments of the type 
worked out in the proof of [SJ Lemma 2]. For details, see Appendix ID. H □ 



4.3 Analysis of Ug, Ug 



4.3.1 Properties of some solution operators 



Lemma 4.3. Let a,a,b G R + . Then 



e~ at (a + bty dt < -e~ ax Y (a + bxf' v [ - X 

i/=0 



oo roo 



x J t 



1 



e~ ar (a + br) 1 drdt < -^^DX! 



v=oe=o 



'i — u)—\ (a + bxY 
a 



Proof. We have, after successive integrations by parts, 
J e~ at (a + btf dt = -e- ax ^{a + bx) 1 



i>=0 

i 



u=0 



i(i-l)---(i-v + l) 
ib 



a 



The statement about the double integral follows from this result. 



□ 



The above lemma can be formulated in the complex plane as follows. 

Lemma 4.4. Let a ,a, b £ M + and v be holomorphic in the half-plane Kez > zq and assume 
\v(z)\ < e - aRe ( 2 ) (a + b\z\y. Then for z withRez > z : 



oo roo 



z Jt 



v(t)dt 
v(r)dTdt 



< -e~ aRe ^y s (a + b\z\) 



3-v 



u=0 
j 3-v 



jb 



a 



v=0l=0 



a 
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Lemma 4.5. Let a > 0, g G R, and / 6e entire satisfying 

l/MI <c f ( q + \z\y^ {z) /orRe(z)<0 

/or some Cf, q, a, a, with (a + a)q > 2j + 1 and < a < a < a. 

XTien the solution u of 

-u"(z) + a 2 u(z) = f(z), u(0) = 5, lim u(z) = 

z— >oo 

satisfies the bound 

, m ^ A/ i 1 i A fe"^ Re(z) /or Re(z) >0 

1 V 71 " ; V« J + 1 |y 7le- oB *W /or Re(z)<0. 



Proof. The proof follows from the appropriate modifications of (6J Lemma 7.3.6]. For details, see 
Appendix ID.2I 

□ 

Lemma 4.6. Let < a < a. Let v be entire and satisfy for some a > 0, C v , b > 0, j G No, the 
bound 

' e -aRc(z) f or Re(z)>0 



\v(z)\ <C v (a + b\z\y 



e -aRc( 2 ) f or Re(z)<0. 



Assuming i— < 1, there holds 



z Jt 



vMdrdt < C, 



1 / 1 



\ aa , 



(a + b\z\y 



' e -aRe(z) / or R e (z)>0 
e -aRe( 2 ) f Qr < q 



Proof. Follows essentially from Lemma 14.41 We sketch the argument for the case Re(.z) < 0. By 
linearity, we may assume C v = 1. We start with the single integral v(t) dt, by selecting as the 
path of integration the line z + r, r G R+, we get with the aid of Lemma 14.41 



v{t)dt 



< g -aRe(z) 



/•—Re z />oo 

/ e-" T (a + 6|z| + 6r) J 'dr + e- aRe(2) / e" aT (a + b\z\ + br) 3 ' dr 

Jo J~Rez 



< e 



-a Ke(z) 



— I 

e^ r (a + b\z\ + br) j dr < e- aRe{z) -^2(a + b\; 



u=0 



< e- SR ^)-(a + b|z|)^( | 6 V < e -*«*)±( a + &|z|)i-* . 

Inspection of the above derivation shows that for Re(z) > 0, the same estimate holds with e - aRe ( 2 ) 
replaced by e _ ^ Re ( 2 ). "We may therefore repeat the argument once more for the function z i— >■ 
v(t) dt to get the claimed estimate. □ 
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Lemma 4.7. Let the entire function v satisfy the hypotheses stated in Lemma\4.b] Then 



\v"(z)\ < 2e a C v (a + b + b\z\y 



-aRe(z) 
-aRe(z) 



for Re(» > 
for Re(z) < 



Vz € C. 



Proof. Follows from Cauchy's formula for derivatives by taking dB\{z) as the contour. □ 
Lemma 4.8. For C\, (3 > and x > the following estimates are valid with 7 = 2max{l, Cf}: 

( Cl t + x) 2e < 2 e (C 1 £) 2i + 2 i x 2e < 1 i (e 2e + x M ), (4.21) 

(4.22) 



supx n e-^ 4x < (— 

x>0 \ e P 



Proof. The result follow from elementary considerations. 



□ 



4.3.2 Regularity of the functions U|l- and 



We turn our attention to equations (|4.15b[) and (|4.15cp . which, after introducing appropriate bound- 
ary conditions, determine ufj, vfj and ufj, vh, respectively. These equations turn out to be systems 
of differential-algebraic equations (DAEs); however, their structure is such that the algebraic side 
constraint of the DAE can be eliminated explicitly and, additionally, we will be able to solve for 4 
scalar functions sequentially instead of having to consider the coupled system. We recall that the 
functions = (ujj,%) T have been defined and studied in Section [IT 



These equations may be solved by induction on j and i. For j = 0, we solve (|4.15bj) for any (i, 0) 
by first solving for uf and inserting it into the equation for vf Q . We have from (|4.15bl 1 st eqn) 



ai 2 (0). 



J i,0 



an(0) - <7j 
which, upon inserted into (|4.15b| 2nd eqn) gives 

f~L\" . Oll(0)o22(0) - O2l(0)oi2(0) „£, 

- [ v i,o) H 



^)E| [4? (0)^,o + 0^(0) 



v i-k,0 



(4.23) 



an(0) 



(4.24a) 



1 " r /o2l(0) aS? ( 0) - a&\o)) ut k , + ( ^||og(0) - a£\o) ) v^_ k 



Ex 
Id 



k=l 



oii(0) 



The above second order differential equation is now posed as an equation in (0, 00) and supplemented 
with the two "boundary" conditions 



wfo(0) = -Ui,o(0), 







as x — > 00. 



(4.24b) 



So, solving (I4.24j) gives us vf Q and then from (I4.23j) we get uf . Inductively, we obtain vf Q and u^ 
for all i > 0. 



11 



Next, we set 



(4.25) 



and we solve with j = (14.15cj, 1 st eqn) for uf (using v\q = 0) with boundary conditions from 



uf o (0) = -Uifi, u^ (x) for x — > oo. 



Then, we solve (|4.15cl 2 eqn) for v^ 2 '- 



Q2l(0)u i (T)dTdt . 



(4.26) 



(4.27) 



In general, assume we have performed the previous steps and we have determined ufj, v^j, ufj, 
f° r au * — an d second index up to j. To obtain the corresponding functions (with j replaced 
by j + 1) we proceed analogously. We first solve (I4.15bl 1 st eqn) for 



u i,j+l 



an(O) ' a n (0) 

and plug it into (|4.15bl 2 nd eqn): 



(4.28) 



" , o 1 i(0)a 22 (0) -O2i(0)oi 2 (0). 



an(0) 



021(0) ,~ £ ,// 
Kj-i) + 



, |Y a 2 i(0) j fc ) 



k=l 



47(o) -o«(o)) 5 



on(O) 



" a n (0) 

021(0) (fc) (&) 



n— fc,j+l 



+ 



on(0) 



a£'(0)-a£(0) 5? 



-fcj+i 



The second order ODE, equation (|4.29ap . is supplemented with the boundary conditions 



(0) = - (t;y+i(0) + S?& +1 (0)) , ^ i+ i(i) -> for 



x — )• 00. 



(4.29a) 



(4.29b) 



This gives us v^ +1 and in turn wfj'+i f rom (|4.28|) . 



Next, we solve (|4.15cl 1 st eqn) for ufj +1 with boundary conditions from Uij+i and ufj +1 : 
- {u^ +1 )" + a u (0)uf J+1 = a 12 (0)v^ +1 - 

min{ij'+l}^ fc 

E 7J («n (0)Sf_ fcJ+1 _ fc - ag (0)«f_ fc J+1 _ fc ) (4.30a) 
fc=l 

-Ki+i(o)+^ i+1 (o)), 

for a? — )• 00. 



^,-+i(0) 



?4j'+i( x ) 



(4.30b) 
(4.30c) 



Finally, we solve (|4.15cl 2 nd eqn) for v^j +3 : 

k {a^(0)uf_ ktj+1 _ k (T) + a^(0)vtk d +i-k(r)} drdt. (4.31) 



v i,j+3\ z ) 



min{i,j+l} 



fc=0 



The following theorem establishes the regularity of the functions ufj, vfj, ufj, vfj. 
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Theorem 4.9. Assume that f , g, and A satisfy \2. 3\) and \2. 7\ ). Let ufj, vfj, ufj, vfj be defined 
recursively as above, i.e., they solve \4.23\) , {^.21$ , {J^.26^ , {1^.25^ , {4-21^ for the case j = and, 
for j > 1 glgp , gfgg j, gWjj , (fCT - Set 

a n (0)a 22 (0)-o 21 (0)012(0)) 

.> u, 
>0, 



a := max | an (0), 
a := min < an(0), 



aii(O) 

011(0)022(0)- 021(0)012(0) 
an(0) 

' e -aRe( 2 ) / or R e (z)>0 
e -aRc( 2 ) j or Re ^) < q 

" L ij' ^ij> U ij> v(< 

Cq, Cy, Ci, Ki, Ki, i = 1, ...,4, independent of e and ^ such that 

^(z)\ < Cl KlK{(Cz(i + j) + \z\) 2(t+3) — 
-' — I tr< f; ;> ~l\2(«+j) 



Exp(z) 



Then the functions ufj, vfj, ufj, vfj are entire functions and there exist positive constants C^, C%, 





»\ 


< 


ft- 


»\ 


< 


Hi 


»\ 


< 




»\ 


< 




= 0. 





1 



Exp(z), 



vt tj (z) < C 2 K* 2 K J 2 (Cv(i + j) + \z\y^>——EMz) 



Exp(z), 
■Exp(z); 



\2(i+j) 



\2(i+j) 



1 

1 

(* + j) 



(4.32) 
(4.33) 
(4.34) 
(4.35) 



Proof. The proof is by induction on j and i. After establishing the claims for the base cases 
(i,j) = (0,0), G {0} x N, (i,j) G N x {0} one shows it by induction on j with induction 

arguments on i as parts of the induction argument in j. The structure of the equations defining 
ufj, vfj, ufj, vfj + 2, is such that one can proceed successively in the induction argument on j by 
providing estimates for v%j, ufj, ufj in turn. In these estimates estimates, one relies on 

Lemma 14.51 for the estimates for vfj, ufj, on Lemma 14.61 for vfj +2 , and on Lemma 14.71 for ufj. 

For details, see Appendix EU □ 



We conclude this section by showing that the boundary layer functions are in fact entire: 

Corollary 4.10. The functions Ufj and Xjfj are entire functions, and there exist constants C , 
7i; 72; P > independent of i, j, n, such that for all x > 

\™th x )\ + \^h x )\ + ISyO*)! + 1^+2^)1 < Ce-^{ +j (i + i)^ 7 2 n Vn G N . 

Proof. Theorem 14.91 already asserted that the boundary layer functions are entire. For the stated 
bound, let n G No, x 6 (0,oo) and use Cauchy's integral theorem for derivatives with contour 
dB n+ \(x). We illustrate the procedure for Uij, the other cases being similar. Theorem 14.91 then 
yields suitable constants C, 7 > such that 

1^)1 ^ C 7^W^7^Tl! iCdi + i) + x + n + l ) 2[l+3) Exp(*)e s(n+1) . 
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With the aid of Lemma \$~8[ we obtain by suitably adjusting C and 7, 

n! 1 



(i + j) 



(c s (» + i) + n + i; 



2(«+.?'),,a(n+i) -az 



Using the observation ((i+j)+n+l) 2(i+i) < (i+j) 2 ( i+ ^(l+(n+l)/(i+j)) 2(i+i) < (i+j) 2(i+i) e 2(n+1) , 
allows us to conclude the proof. □ 



Corollary 14.101 shows that the terms defining the boundary layer contributions and TJg L are 
indeed of boundary layer type. A summation argument then shows that also U|f L and have 
this property provided Mi and M2 are not "too large", viz., M\ = 0(fi~ 1 ) and M2 = 0((^/e) _1 ): 

Theorem 4.11. There exist constants C, 5, 7, K > 0, independent of e and \i, such that under 
the assumptions n{M\ + 1)K < 1 and Ej\i{Mi + 1)K < 1, there holds for the boundary layer 
functions V^l an d ^bl °f B3P an< ^ EZP> ^at, upon viewing U^f L and \J^ L as functions of 
x (via the changes of variables x = x/fi, x = x/e etc.), we have G BC°°(5^, C, 7) and 

Ufii £ BC aD (5e, C, 7). Furthermore, the second component v of\Jg L satisfies the stronger assertion 
v G^£°°(fe,C(e/^) 2 ,7). 



Proof. We do not work out the details here since structurally similar arguments can be found, for 
example, in the proofs of p3 Thm. 3] or JpJ Thms. 7.2.2, 7.3.3]. Essentially, by inserting the bounds 
of Corollary 14.101 in the sums defining U|^, Ug L and using the conditions \i{M\ + 1)K < 1 and 
ej \x{M2 + 1)K < 1 for K sufficiently large, one obtains upper estimates in the form of (convergent) 
geometric series. We point out that the sharper estimates for the second component v of Ug L 
stems from the fact that V{ = v% \ = 0. □ 



4.4 Remainder estimates 

In this section, we analyze Rm- This is done by estimating the residual L £j/j Rjv/ and then appealing 
to the stability estimate (12. 9p . We will estimate L £jA1 Wjv/ — F, L £jfl XJg L , and L £iA1 U^f L . 

4.4.1 Remainder resulting from the outer expansion: control of L £ „Wm — F 



Theorem 4.12. Let U be the solution to the problem (12, ip . Then there exist 7, C > depending 
only on f, g, and A such that the following is true: If M\, M2 G N are such that \xM\^ < 1, then 
with Wm given by (j4.3j) we have 



IL^ (U-W M )|| Loo(/) <C7z 2 



H-e V/V 



\ m 2 +2 
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Proof. We have 

(U - W M ) 



Mi M 2 , ,j 
Mi M 2 

i=0 j=0 vrv 



-e 2 ^- + annjj + ai2%' 
-^ 2 ^- + a 2 iUij + a 22 % 



Defining the sets 



= : * < Mi, j < M 2 , i > Mi - 1 V i > M 2 - 1}, 

/« := {(i,i) : i < Mi, j < M 2 , Mi - 1 < i < M x }, 

we see, after some calculations, that (|4.15ap and Lemma 14.21 imply 



(4.36) 
(4.37) 



L e „ (U - W 



My 



Mi M 2 



i=0 j=0 1 



'J 



5ij 



+ 



a i-2,j-2 £ u ij 



IJ 



hi 



Hence, with the aid of Lemma [4.2l and Cauchy's Integral Theorem for derivatives, we get for a fixed 
5 > in the statement of Lemma 14.21 



\L E!fl (U-W M )\\ L oom < CC S 



Mi -2 



M 2 



Mi 



M 2 



^Eqm- 1 *)* e ( £ /M) i+2 +M 2 E (M-^rE^/^ 2 

i=0 j=M 2 -l i=Mi-l j=0 

Mi M 2 

+/" 2 E (^^T^) 3 ' • 



i=Mi-l 

Hence, by selecting 7 = 5~ 1 K/2 we get 



i=o 



|X £)M (U - Wm)IUoc (j) < CC S ^ - + (MM17) 



M 2 +l 



\Mi-l 



1 



l-e//i 



£ 



M 2 + l 



+ (^M l7 ) Ml 



Mi 1 



[i — e 



□ 



4.4.2 Remainder resulting from the inner expansion on the e-scale: L £ifl Ug L 

We next consider the inner expansions. We will only consider the contribution from the left 
endpoint as the contribution Vj^ from the right endpoint is treated completely analogously. To 
simplify the notation, we drop the superscript L in ufj, vfy 

In order to simplify the ensuing calculations, we employ the convention that 

Uij = Vij = for i > Mi or j > M 2 and A k = Vk < 0, (4.38) 
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and let the summation on i and j in the definition of run from to oo. We recall that the 
differential operator L £jfJi takes the form (|4.10p when applied to functions depending solely on x, 
and compute L £ JUbl ( c ^- 



for the definition of V^ L ) 



T TT M 



i=0 j=0 




„ + EEEf* 



i,i+2 



i>0 j>0 k>0 



+ ^ X k A h 

k>0 



£V(£ 



Vi—k,j—k 




H k x k A k 



where the fact that u^o = Ufc,! = was used. We see that "equating like powers of fi and e//x " 
yields equation (|4.15cp . hence there will be no contribution to the sums for when both i = 0, . . . , Mi 
and j = 0, . . . , M2 — 2. Moreover, the convention (|4.38p implies the following restrictions on the 
sums: 



k < min{i, j}, 



fc<Mi, j-k<M 2 . 



(4.39) 



Thus, we have 



r frM 



E 



ij':i>Mi+l 
or J>M2 — 1 



AV 



min{ij'} 

E 

fc=max{i— Mi J—M2} 



x k Ai 



V-i—kJ—k 
Vi—kJ—k 



Mi M 2 

EE"' 

i=0 j=M 2 -l 



(4.40) 







Using the estimates (14.34p . (14.35P for Uij, Vij of Theorem 14, 9} one can estimate this triple sum to 
obtain the following result for the remainder on the positive real line: 

Theorem 4.13. There exist C, 7, K > depending only on A, /, g such that under the assump- 
tions 

< jxe < 1 and /x(Mi + 1)7 < 1 and -(M 2 + 1)7<1, 

A* 



we /iaue 



< Ce'^l 2 {[K{M X + l)/x) Ml + (K(M 2 + lje/zx)^" 1 ) . 



For 7x6 > 1 we have, under the same conditions on (Mi + l)/x and (M2 + \)ej \x, that 



< Ce~^ /2 . 



Proof. The proof relies on using the estimates (|4.34p . (|4.35p that are available for and V{j. The 
triple sum in (|4.40p is estimated by using convexity of the function k 1— > j k (i — k) l ~ k (j — k)i~ k and 
k 1 — y k k j k (i — k) l ~ k (j — k)i~ k and by considering the two following two cases separately: 

(i > Mi + 1 V j > M 2 + 1) A (i - Mi < j - M 2 ) A (j - M 2 < k < i) A (i < j), 
(i > M% + 1 V j > M 2 + 1) A (i - Mi < j - M 2 ) A (j - M 2 < k < j) A (j < i). 

For details, see Appendix ID. 41 □ 
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4.4.3 Remainder resulting from the inner expansion on the ^-scale: L e uUbl 

In a similar fashion we may treat the left inner (boundary layer) expansion associated with x = x/fi 

(cf. 631)), 

where we have dropped the superscript L for notational convenience. We recall from fj4.9j) that 
for functions U depending only on the variable x, the differential operator L £tfl takes the form 
L £i/ jXJ = — /i _2 E £,M U" + Yl'kLo I 1 % AfcU. In order to simplify the ensuing calculations, we employ 
the convention that 

Uij = Vij = for i > Mi or j > M 2 and A k = VA; < 0, (4.42) 
and let the summation in (|4.4ip run from to oo for both i and j. We calculate 

co oo /II \2~tl \ 

A—c\ a—c\ \ 13 ' »>>r\ u^n 



i=0 j=0 y 7 i>0j>0k>0 

-u'j 

" %3 7 fc>0 



i>0 j>0 



- V ij J t?n V Vi-kj 



where the convention Uk-2 = Uk,-i = was used. As expected from the derivation of (|4.15b|) . 
"equating like powers of fx and s/fJt" yields eqn. (|4,15bp . hence there will be no contribution to the 
sums for when both i = 0, . . . , M\ and j = 0, . . . , M 2 . Moreover, the convention (I4.42p implies the 
following restrictions on the sums: 

i < M\ for the terms involving u'-j_ 2 , v"j, (4.43) 

< i - k < Mi for the sum on k, (4.44) 

j < Mi for the terms involving u%j, v iji u i—k,j> v i—k,ji 

(4.45) 

j < M2 — 2 for the terms involving u'- j_ 2 - (4.46) 
Hence, we arrive at 

Ls,Ml = (4-47) 

00 M2 / \ j i / ~ \ Mi A/2+2 , \ j / ~ii 

e £"■(£) e +E e { - u r 

i=Mi+l j=0 VP/ k=i-Mx V 7 i=0 j=M 2 +l VP/ V 

Using the bounds of Theorem 15.51 these sums will be estimated, when x > 0, in the following: 
Theorem 4.14. There exists C, 7, K > smc/i i/mf under the assumption 

< < 1, and 7Ai(M x + 1) < 1 and ~7(M 2 + 1) < 1, 

we aaue 

L e , M U^(x)| < Ce"^/ 2 ((K(M X + %) M i +1 + (tf(M 2 + l)e/^ +l 
For Jlx'y > 1 and the same assumptions on fi(Mi + 1) and ej n{Mi + 1), we have 
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Proof. The proof is split into two cases: for xfij > 1, one exploits the fact that Ug L (and its 
derivatives) is exponentially small. For the converse case x/i'j < 1, one bounds the sums (|4.47|) . 
For details, see Appendix ID.5I 



4.5 Boundary mismatch of the expansion 



□ 



Theorem 4.15. There exist constants C , b, 7 > such that under the assumptions 

e 



H{M X + 1)7 < 1 



and 



-(M 2 + l)7< 1, 



one has 



\w M + uf L + u^ L + v# L + vMlWl-od < c 



e -b/e + e -b/n 



Proof. For Rm(0) at the left endpoint, we note that U(0) = gives 

R M (0) = U(0) - \W M (0) + Uf L (0) + Uf L (0) + V&(l/ M ) + Vf L (l/e) 



The condition (I4.15dj) for the boundary conditions of the left endpoint for the terms Uj j and Ujj, 
produces W M (0) + Uf L (0) + Uf L (0) = 0. Hence, it remains to estimate |Vf L (l//i)| + |Vjf £ (l/e)| 
which can be done based on Theorem 14,91 □ 



4.6 Proof of Theorem 14.11 



From the estimates for the residual L £j/j Rm of Theorems 14.121 14.131 14.141 we infer the existence of 
q > such that under the assumption 



fi < q and 
the choice Mi ~ 1/fi and M 2 ~ fi/e yields 

\\L £ ,hRm\\l°°(i) < C 



e 

- < 



-bfi/e 



where C, 6 > are independent of and e. Hence, by stability and Theorem 14.151 we get that 
IIR-mIIe,/ < C [e^ b ^ 1 + e~ 6 ^/ e l . The sharper result for the second component u of U follows from 
the fact that v i:0 = = 0. 

It remains to formulate a decomposition of U for the case that fj,>qoxe/fj,>q. In this case, we 
have e~ h ^+e- b ^ /e is 0(1). Given that \\U\\ E ,i = 0(1), the trivial decomposition U = 0+0+0+R M 
provides the desired splitting. 
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5 The first two scale case: Case IIIII 



In this case it is assumed that is not deemed small but e/fi is deemed small. The main result 
of this section is the following regularity assertion: 

Theorem 5.1. There exist constants Cyy , ^yw, Cbl, Ibl, $, b, q > independent of e and /j, 
such that for e/fi < q the solution U of h2.1\) can be written as U = W + TJbl + R> where 
W € A(fi,Cw,1w), Url G SC°°(5£,Cbl,Jbl)- Furthermore, R satisfies 

+ ll^e,M^'lli 00 (-/') — C e ■ 

In particular, ||R||£,/ < Ce~ h l e . 

Additionally, the second component v of TJ bl satisfies the sharper regularity assertion 

5eBr(fc,C Bi (^f )7BL ). 

We employ again the notation of the outset of Section H] concerning the stretched variables x = x/e 
and x = (1 — x)/e and make the formal ansatz 

oo 



/VV r - i 

U (^)~E - U,(x)+Uf(x) + Uf(x fi ) . (5.1 



We proceed as in Section|4]by inserting the ansatz (15. lh in the differential equation (|2.1ap . separating 
the slow (x) and the fast (x and x ) variables and equating like powers of e/fj,. We also recall that 
the differential operator L e ^ takes the form (|4.10p on the x-scale. The separation of the slow and 
the fast variables leads to the following equations: 

(-E^U? + A(s)U<) = F =(^)> ( 5 - 2a ) 

W£ -e-^Ufr + 5> fc - x fc A fc Uf = 0, (5.2b) 
•=o w V fc=o w / 

and an analogous equation for U^. Writing again 

tWM, of=(S 



we obtain from (|5.2p by equating like powers of e//x: 



(5.3a) 
(5.3b) 

(5.3d) 

k=0 



k=0 



- n 2 u'-_ 2 + a u ui + a 12 Vi 


— fii 




= 9i, 


4?(0) si , ag>(0) L \ 

fc , «,-fc 1 k] V i—k j 


= o, 


<&<p)*l , 4?(o).l A 

it! 1 kl Vl ~ k j 


= o, 



19 



where we employed the definition of A&, the notation /q = /, g$ = g as well as /j = g$ = 
for i > 0, and the convention that function with negative subscripts are zero. Corresponding 
equations are satisfied the functions uf and vf. The boundary condition (|2.1bj) is accounted for by 
stipulating Uj(0) + Uf (0) = and U»(l) + Uf (0) = for all i > and suitable decay conditions 
for \jf as x — > oo and, correspondingly, for as x R — > oo. Rearranging the above equations and 
incorporating these boundary conditions, we get a a recursion of systems of DAEs in which the 
algebraic constraints can be accounted for explicitly. We obtain for i = 0, 1, 2,... : 



-mV + (a22ail a ir a2l) ^ = ft - tr (/< + A* 2 <- a ) 

= -£f(0), u i (l) = -u*(0) 



«i(0) 



an 



: )" + aii(0)uf 



k=l \ 

u\ (0) = -Uj(0) , uf -)• as x ->■ oo 



aia(0)^ 



'4?(o) 



with 



k=0 



A;! 



= Si 



, «_ 



/i 



/ if i = 0, 

otherwise ' ' 



i = V i > 0, 

o if i = 0, 
otherwise 



(The functions u R ,v R satisfy similar problems as (|5.6p and (|5.7p . respectively.) 



(5.4) 
(5.5) 
(5.6) 

(5.7) 



5.1 Analysis of the functions Uj, Uf , U 



R 



5.1.1 Properties of some solution operators 



satisfies 



lk {n) ||L-(/) < Cfmaxln,^ 1 }" [C 5 + \g+\ + |o_|] Vn G N , 
/or a constant C that depends solely on the function c. 



(5.8) 



Lemma 5.2. Assume that the function c is analytic on I and c > c > V z £ 1. Let ^ G (0, 1]. 
Then there exists 70 > (independent of fx) such that for all 7 > 70 and all C g > the following 
is true: If g satisfies 

||2 (n) l|L~(/) < a^maxK^ 1 }" Vn G N , 
then the solution u of the boundary value problem 

—fj?u" + cu = g in I 
u(0) = o_ G E , u(l) = o+ G M 



20 



Proof. This follows by inspection of the proof of [5J Thm. 1] if one replaces the energy type argu- 
ments by an appeal to the comparison principle to get L°°-estimates. 

For details, see Appendix lB.il □ 

Lemma 5.3. Let m E No and g be a function analytic on I. Then there exist C , 70 > depending 
only on g such that the following is true: If v satisfies 

\\v {n) \\ L oo {I) < C V H max{n, /i" 1 }" Vn € N , 

for some C v , fi > and 7^ > 70, then the function V := gv^ satisfies 

||^ (n) || L -(/) < C'a7™ +m max{n + m,/i- 1 } n+m Vn E N . 

Proof. Since g is analytic, there exist C g , ^ g > such that 

<C gl n g n\ Vn G N. 
By Leibniz' rule, we get for 7^ > 70 > 7 9 e in view of v v < v\e v : 



< C v C g il V" - n ' ( ) max{n + m, n~ 1 } 



—l\n+m—v 



1 - e7 9 /7„ 



Lemma 5.4. For < S < l/(2e) and i 6 N there holds 

f-i±* r < 2, 



5.1.2 Regularity of the functions Uj, Uf, Uf 

Theorem 5.5. Lei Vi,Ui,u^ ,v\ , uf, be the solutions of ()5.4j) (|5.7[) . respectively. Let G C C be 
a complex neighborhood of I and set (3 = y^O]T(0) G ^- T/ien i/ie functions uf and v\ , uf- and vf- 
are entire, and the functions Ui, V{ are analytic in a (fixed) neighborhood of I. Furthermore, there 
exist positive constants 7, C u ,C v ,Cp, Kp, £ = 1, 4 independent of e and [i, such that 

= v\ = 0, (5.9) 
\vf +2 (z)\ < C X K\ (u + JL^ ~ (iC v + \z\ fe-^ z \ (5.10) 



□ 



Proof. See Appendix IB.21 □ 
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(z)\ < C 2 Ki L + i (»C„ + M) 2i e -W*), (5.11) 



- (n)|l L-(/) < ^((i + nj/i + lfVmaxjn./i- 1 }" Vn G No, Vi G No, (5.12) 



< C 4 ^i((i + n)/i + l) i ^7 n max{n,^ 1 } n Vn G No, Vz G No- (5.13) 



Furthermore, analogous estimates hold for uf and with f3 being now y/an(l). 

Proof. The proof is by induction on i. The general structure is to estimate Vi, Ui, Ui, and Vi+2 in 
turn using Lemmas 14.51 14-61 15.21 15.31 For details, see Appendix IB. 31 □ 

We conclude this section by showing that the boundary layer functions are in fact entire: 

Lemma 5.6. Let (3 > be as in Theorem 15.51 (i.e., (3 = \J an(0) ). The functions XJf are entire 
functions and there exist constants C, 71, 72 > independent of i, j, n such that for all x > 

\u^\x)\ + \v$ 2 (x)\ < Ce-M^Wn + 1)' 7 2 Vn G N - 
Proof. Similar to the proof of Corollary 14.101 □ 
By a simple summation argument we get from Lemma 15.61 



Theorem 5.7. Let (3 be as in Theorem \5.5\ (i.e., (3 = y/an(0) )■ There exist constants C , 7, K > 
such that for 7{(M + l)e + — } < 1 the function 



M / \i 



8=0 

satisfies for x > 



U 



M ^ W 
BL 



) n (x)\< Ce-P /2x K n Vn G N - 



^4n analogous result holds for Vj^ L := X^f=o( e /^) i ^? ™^ fi^fen by f3 = y/au(l). 

Proof. See Appendix lB~4l for details. □ 

5.2 Remainder estimates 



We now turn to estimating the remainder obtained by truncating the formal expansion (15. ip . We 
write 

U(X) = ( "(x) ) = Wm{x) + ^ (5?) + ^(2) + R ^)> ( 5 - 14 ) 
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where 

M 



denotes the outer (smooth) expansion, 

denote the left and right inner (boundary layer) expansions, respectively, and 



R M (x) := ^ j = U(x) - (W M (x) + U^(x) + V£ £ (s«)J (5.17) 

denotes the remainder. Our goal is to show that the remainder Rjv/(x) is exponentially (in e) 
small. First, we need to obtain results on the other terms in (|5.14p . Note that by the linearity of 
the operator L £ ^ we have 

L £ ^R M = L e ^ (U - W M ) - L £ ^ L - L £ ^ L . 
The terms on the right-hand side are treated separately. 



5.2.1 Remainder resulting from the outer expansion: F — L f)/1 Wj,f 
We have the following theorem. 

Theorem 5.8. Let U be the solution to the problem (|2.ip . Then there exist C , 7 > independent 
of fx, e, and M such that for given by (|5,15[) we have 



\L £ ,, (U-W M )|| ioo(J) <C 



_ \ M+l 

7^) +( 7 (M + 1) £ ) M+1 



Proof. Using (|5.2a|) (j5.4|) . ()5.5[) we obtain, after some calculations, 

||^(U-W M )|| LOO(/) = £ (- [i 2 ( ^ 



i=M-l 



From Theorem 15.51 we therefore get, with the observation ^(M[i + 1) = Me + ^, the desired 
result. □ 



5.2.2 Remainder resulting from the inner expansion: L £ ,^U^f L 

We now turn our attention to estimating L £r/J TJg L . Since L e/t V^ is treated with similar argu- 
ments, we will not work out the details. We have: 
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Lemma 5.9. The functions U|fr satisfy 



t fr M 



E 

i>M+l 



i i 



M 



E V k x k A k XJt k + e 



k=i-M 



i=M- 



1 V fe)" 



(5.18) 



Proof. In order to simplify the calculation of L £ ^XJg L , we employ temporarily the convention that 



uf = vf = Vi > M + 1. 



(5.19) 



With this convention, we calculate (cf. (jOcj) . (JSHd])) 



r tt m 



E 



i=Q 



/i 



fc! n i-fc "i fc!^ 



V 



fc] "i-fc fc! y i-fc 



/ 



£)E^A fc uf_ t - E (j)E 

p/ fc=o i=M-l VP/ fc=0 



4^(0)^ ; 



i>M+l vr " / fc=0 i=M-l sr ~ / fc=0 \ fc! i — k k\ u i-k 

here, we employed additionally (|5.3c|) . (|5.3dj) to see that the terms corresponding to i G {0, . . . , M 

< % - k < M 

so that we obtain 



2} in the first sum are zero. Finally, our convention (|5.19p and the fact that = = for j < 
implies the restrictions 



£) Y,^ k ^tk= E m E ^A*©f 

.. .w-i-i ^' fc=0 i>M+l k=i-M 



-k ■ 



which produces the double sum in (15. 18H . Lifting now the convention (|5.19p we can use (j5.3c[) to 
replace £* =0 ^ k (^uf_ k + ^T^f-fc) with -vf +2 {x). □ 



It remains to bound L e ^\Jjf L . For that purpose, we need the following lemma: 

Lemma 5.10. Assume < xe^a < X < 1 with ja, X, C\, (3 > known constants. Let < e < 
fj, < 1. T/ien 



1 V fc 1 



(Ci(i-it) + x) 2 ^e-^ 



M+l -/3x/2 



< r ^(7£(M + l + l/ M )) JM+1 e 
where C, 7 > are positive constants independent of e and \i. 



Proof. The key ingredients of the proof are the estimates given in Lemma 14.81 For details, see 
Appendix IB. 51 □ 
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Theorem 5.11. Let (3 > be given by Theorem 15,51 i.e., (3 = y/au(0) There exist constants C, 
7 > such that, under the assumptions 



e{M + l) + - ) 7 < 1 £ J . 



U given by f|5. 16j) satisfies 



Le,^ L {x) < C(s(M + 1 + l/vhr-'e 



For the case xe^f > 1 (but still assuming \e{M + 1) + - ) 7 < 1), we have 



< Ce'P/ 2 *. 



Proof. For the case of xe being sufficiently small, the starting point is (|5.18p . which is a sum of two 
contributions, namely, the double sum and the single sum (consisting of merely two terms). For 
the double sum, we recall that x = x/e, so that Lemma 15.101 produces the desired estimate. For 
the single sum, we use the estimates of Theorem 15.51 for vh 2 . From Cauchy's Integral Formula for 
derivatives with contour being taken as dB\(x), we get 



\{v, 



(i + 2) 



■(C(i + 3) + \x\) 



Using (I4.21D and (I4.22D . we see that this contribution can be bounded by C*( 7 (M + lje/n)* 1 - 1 for 
suitable C, 7 > 0. 



For the case xej > 1, we use the exponential decay of expressed in Theorem 15.71 The factors 
e~ 2 , which arise when computing L £ „Um can be absorbed by the expontially decaying term since 
< 7X (see Lemma l4.8p . □ 



5.3 Boundary mismatch of the expansion 



Theorem 5.12. There exist constants C , b, 7 > such that under the assumptions 



one has 



e(M + 1) + - 7 < 1, 



||W M + + VM L \\Lo° m < Ce-V £ . 



Proof. We consider the left endpoint of I (the right endpoint is similar). By construction, 



|Rm(0)|| 



U(0) - (W A/ (0) + U^(0)+V^(l/ £ 



rM 



V M (l/e) 



The result now follows from Theorem 15.5 



(5.20) 
□ 
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5.4 Proof of Theorem 15.11 



The result of Theorem 15.11 now follows from combining Theorems 15.81 15.114 15.121 



6 The second two scale case: Case IIVI 



Recall that this occurs when /i/l is small but e/fj, is not small. The main result is: 



Theorem 6.1. The solution U of \2.1\) can be written as U = W + XJbl + R; where W € 
A(l,Cw,Jw) and ^Jbl £ BC 2 (5/j,,Cbl,1bl), for suitable constants Cw , Cbl, Iw, Ibl, 5 > 
independent of [i and e. Furthermore, R satisfies 

||R|| Loo(a/) + ||L £iM R|| L2(/) < C(n/e)- 2 e~ b ^, 
for some constants C, b > independent of fi and e. In particular, ||R||£j < (/i/e) 2 e _f '/ M . 



In this case we recall the stretched variables x = xj\x and x = (1 — x)/[i and make the formal 
ansatz 

oo 

u 

... (I 

where again 



oo 

J>< [UiW + Of (J) + uf(J")] . (6.1 



U iW .(«), Uf ® - ( ||> , . 

and an analogous definition for XJ R . We also recall that the differential operator L £jIM takes the 
form (|4.9p on the x-scale. The separation of the slow (x) and the fast variables (x and x R ) leads 
to the following equations: 

E ^ (-E £ '"U? + A^Ui) = F=(0, (6.2a) 

oo / oo \ 

-/i- 2 E^(Uf)" + E^A fc Uf = 0, (6.2b) 



i=0 \ fc=0 
jR 

powers of [i to get with the matrix 

2 



and an analogous system for U^. Next, for e appearing in (|6.2p we write e = \i~ and equate like 







5_ 

I 



the following two recurrence relations: 

-E^Uta + A^U* = Fi, (6.3a) 

-E*/* 1 ^)'' + A*_ fc u£ = 0, (6.3b) 



k=0 
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where, as usual functions with negative index are assumed to be zero and Fj are defined by Fq = 
(f,g) T and Fj = for i > 0. The terms of the outer expansion, Uj are obtained immediately from 

- CM^Y 

The functions Tjf of the inner expansion are defined as the solutions of the following boundary 
value problems: 

- E^Uf )" + A(0)Uf = -J2^~ n Al n Vt (6.5a) 

n=0 

Uf (0) = -Ui(0), Uf(x) ^ as x — > oo, (6.5b) 
and an analogous system for U^. 



6.1 Properties of some solution operators 

The functions XJf of the inner expansion are solutions to elliptic systems. In contrast to the 
previous arguments, which were based on estimates for scalar problems (for which strong tools 
such as maximum principles are readily available), we employ more general energy type arguments 
here to deal with the case of systems. We start by introducing exponentially weighted spaces on 
the half-line (0, oo), by defining the norm 



oc 



|| u ||2 := / e W*\ u {x)\ 2 dx, (6.6) 

Jx=0 

with the obvious interpretation in case u is vector valued. The following lemma shows that elliptic 
systems of the relevant type (|6.5p can be solved in a setting of exponentially weighted spaces: 

Lemma 6.2. Let v G (0,1] and set E := f ^ 1 1' B £ ^ 2 positive definite, i.e., 

x T Bx > /3g|x| 2 for all x G M 2 . Then the bilinear form 

/•oo 

o(U,V)= / U'-EV' + U-BVds, 

Jx=0 

satisfies for a constant C > that depends solely on f3$, the following inf-sup condition for all 
< j3 < j3 : For any U G ^(0,oo), where 

i^(0,oo) = {u : ||u||i >/3 < oo}, (6.7) 
there exists V G H-b w ith V/0 such that 

V) > C__^_[|XJ]| 1 ^]|V]| lj _ ;S . 

Here, we define for a G K 

/■oo 

||U||? := / e 2ax [U' • EU' + U • BU] dx. (6.8) 

Jx=0 
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Proof. Given U G HUO,oo) we set V := e 2 ^ x U. Then V'(x) = 2(3e 2 P x XJ(x) + e 2 ^ x V'{x) and thus 



l|V||U < 4(1 + M) 2 ||U||^ 
o(U,V) = ||U||?fl + 2P r e 2 ^U'-EUdx> ( 1 - ^ ]|XJ]|f ^. 

The result then follows. □ 

Lemma 6.3. Let f satisfy ||f||o /? < oo for some (3 G [0, /?o), afid fe£ g G M 2 . Then the solution U 
of 

-EU" + BU = f , U(0) = g, 

satisfies for a C > independent of (3, the estimate 

[IXJIU,^ < C7(A) - []|f [|o,>s + |g|] • 

Proof. Let /3i > /3q and set Uo = ge~@ lX . Then Uo satisfies the desired estimates. The remainder 
U — Uo satisfies an inhomogeneous differential equation with homogeneous boundary conditions at 
x = 0. Hence, Lemma 16.21 is applicable and yields the desired result. □ 

Lemma 6.4. There exist 5q > and Co > such that for every 5 € (0, 5q] and every m G No, 
there holds 



{i + mf (2n + 1 + m) 2n+1+m ~ °' 
Proof. See Appendix lC.il □ 

6.2 Regularity of the functions Uj, Uf , Uf 



Lemma 6.5. The function Uj defined by fj6.4\ ), are holomorphic in a neighborhood G of I and 
satisfy for some C, K > 0, 

\Ui(z)\ < CS^KV VzeG 5 :={zeG\ dist(>, dG) > 5}, Vi > 0. 

Additionally, U2«+i = /or a// i G No- 

Proof. We note that < 1. The arguments are then analogous to those of Lemma 2]. The 
arguments are also structurally similar to the more complicated case studied in Lemma 14.21 □ 

We now turn to estimates for the inner expansion functions Vf . 

Theorem 6.6. The functions Uf defined by \6. 5a\) . \6. 5b\) are entire functions and satisfy for all 
(3 G (0, /3q) (with (3q given as in Lemma W^ . 

||Uf ||i i/3 < CK\Pq - /?)-( 2m ^ Vi G N . (6.9) 
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Proof. The case i = follows from Lemma [6,3l and Lemma [6, 51 For % > 1, we proceed by induction. 
In order to be able to apply Lemma 16.31 we define 



2-1 



F(x) = Y J ^i-nX i - n U L n (x). 



n=0 



Next, we estimate, for an arbitrary ft E (0, /3q) and (3 E (/?, /3q), with the aid of Lemma l4~ 

/»oo 

e^^-^lU^S)! 2 ^ < / e-W-^x 2 ^- n h 2 ^\ijL(x)\ 2 dx 



Jo 

< supe-^-^x^HU^II^ 



< e -2(i-n) 



i — n 



P- 



2(i-n) 



|U 



Li|2 



'0,(8 



< C 2 ^ 2n n 2 ™(/3 " ^)-2(2n+l) g -2(i-n) 

Selecting /? = (/?o — + /3 for some k G (0, 1) to be chosen shortly, we get 

e 2 ^x 2 ^\V^(x)\ 2 dx< 

1 



i — n 



P-/3 



2(i-n) 



< C 2 K 2n n 2n (f3 - / 3)-2(i+n+l) e -2(2~n)(. _ „)2(t-n) 



K 2(i-n)(i - K ) 2 (2n+1) 



The choice k 



i+n+l 



yields 



/ e 2/3 ^ 2(i - n) |U^(x)| 2 dx < C 2 K 2n ((3 - /3)- 2 - r 
Hence, 

i-l 



(n J- 4 J- ~\ \2(i+n+l) 
{ i+n+l )f ,-2(i-n) (* + I + 1) ^ } 

(2n + 1 2 (2™+ 1 ) v ; 



n=0 



i» (2n + l) 2 «+i 



n n (i + n+1) i+n+l 

j 'A - " 7KZ I T^Ti (PO-Pj 

) 

< CC^Co^o-/?)- 2 '^- 1 ^, (6.11) 



where we appealed to Lemma [63] and used implicitly that i<C is sufficiently large. Using Lemma [63 
for a fixed 5, we get from Lemma [67 



|U;[| li/3 < (Po-P)- 1 [\\F\\ 0jP + C(K/Syi 
< CiPo-p)-*- 1 !?? 



The expression in square brackets can be bounded by 1 uniformly in i and /3 6 (0, /3q) if we assume 
that C and are sufficiently large. □ 



We next refine the argument to include bounds on all derivatives of Uj : 
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Theorem 6.7. There exist C\j , K\, K2 > 0, independent of (3, v, m and i, such that 
II (Uf ) (m) ||o l/3 < Cutfa - P)-^ +l+m \i + mfK\K^v- m . 

Proof. The cases m = and m = 1 are covered by the above lemma (note: ||E _1 || < v~ 2 ). The 
remaining cases are obtained as usual by differentiating the equation satisfied by Uf and then 
proceed by induction on m. For details, see Appendix EU □ 

We conclude this section by showing that the boundary layer functions 

M 



1 L 

'1 1 

8=0 



are entire functions. 



Theorem 6.8. Fix f3 € (0, (3$). There exist constants C , 7, K > such that under the assumption 
H(M + 1)7 < 1 there holds 



dx' r 



fjM 



< CK m u- m Vm G N . 



Proof. From Theorem 16.71 we see that for all m 6 No, 



M 



ii(u^) (m) iio,^ < c(A-r 1 E^ ! ('+mr^-/?)- 2 w^ m 

i=0 

M 

< C{fio - $)- x K%v- m £(2(/3 - py 2 Kx^y + (2Kinm,y 

i=0 

< C(/3o-/3)- 1 iv->- m ^(2K 1 ^M)M(/3 -/3)- 2j +g) j < Ctf 2 m 



i=0 

for an appropriate K2 (depending on (3q — /?!), if we assume that fiM is sufficiently small. The key 
observation for this fact is to note that for m > M we have 

M M -I- 1 M 4- 1 / m \ (M-l)A")m 

£(m/M)* < (M|l)WM) J, = m ^WMr<m^ 

i=0 ^ ' 

and n 1 /™ — > 1 for n — > 00. □ 



6.3 Remainder estimates 

6.3.1 Remainder estimates for the outer expansion: F - I £i/i W^ 



As before, the formal expansion (16. 1|) is truncated after M terms to yield the decomposition 

M M M 

U = XyUiCs) + £//Uf (5?) + £//Uf (S fl ) + R^ 



,M 

fj, ) T if 

i=0 i=0 i=0 
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A calculation shows 



M 



F - L £iM E //U, = ^ M+2 E^UX f . 



i=0 



We therefore get 

Theorem 6.9. There exists 7 > independent of e, [i such that for fj,(M + 1)7 < 1, there holds 

M 



|F - L £ ^ J2 |k°o ( ,) < C(/i(M + 1)7) 



M+2 



8=0 



Proof. The proof is analogous to that of Thm. 6] . 



□ 



6.3.2 Remainder estimates for the inner expansion: L £ ^XJg L 



We consider only the contribution from the left endpoint and define := Yli^o f^^f ■ A 
calculation shows 

M 

Le,M ! L = E ^E^ fcA *-^- (6-i2) 

i>A/+l fc=0 

The following lemma provides an estimate (in an exponentially weighted space) for L £ ^\J^ L near 
the left endpoint: 



Lemma 6.10. There exist C , 5, f3 > 0, K > such that 

2 

-- x - V" „' V" ;.'-/•• a , tt l 1 d ~ < C (Kn(M + 1)) 2 ( M+1 ). 



5/V 



x=0 



E m'E^a^u^s 

i>M+l fc=0 



Proof. For fixed x > we estimate 

M 

E ^E ?_fcA -^(5 

i>M+l fe=0 

For ^X7 A < 1/2, we estimate further 



M 



<C A E |U£(x)| E 7^*5*- V 



fc=0 



j=M+l 



3.13) 



E 7a *2 i ~V < 2(/ii7 A ) M+1 (i7A)" fc < 2^ M+1 (i 7A ) A/+1 - fc . 
=M+1 



Inserting this in (|6.13p . we see that employing the estimate (|6.10p we can reason in exactly the 
same way as we have to reach (|6,lip . to get 



5/11 



x=0 



M 



;>M+1 fc=0 



dx < C(Kfi(M + 1)) 



2(M+1) 



which is the desired estimate. 



□ 
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Lemma |6. 101 provides an estimate for L e „U^ near the left endpoint; the following result provides 
an estimate on the whole interval I: 

Theorem 6.11. There exist C , 7, K , b > such that for n(M + 1)7 < 1, there holds 



(/i(M + 1)K) M+L + 



M+l + e -b/fi 



Proof. We merely consider the left boundary layer. Lemma [6.101 allows us to estimate for L £ /jXJMt 



on the interval (0, 5) with the change of variables x = \xx: 

Jx=0 



[ \L £ ^ L \ 2 dx = l x\ \L £ ^XJ A J L (x)\ 2 dx < Cfi((M + l)fiK) 

Jx=0 



2(M+1) 



For the interval (5,1), we note that 



L e>lt V% L = — E^' 1 ( TJ# L ) (x) + A(x)U% L (x). 



\L e ^ L \ 2 dx < 



Hence, we can estimate 

J x=S Jx=6/[i 

With Theorem 16. 8\ we therefore arrive at 

T M 1 1 2 



C/m / I (uf L ) (x)\ 2 +\V^ L (x)\ 2 dx < Ce~ 2 ^ 

Jx=S/u V 7 



Tp IP 1 I TTM ip 

U BL I II 0,0 + II u BLIIo, / 



2(M+1) , „-4--2*j8/a» 



II^UfLlll 2 (/) < C^^(M + l)K)^^+^e 

□ 

Remark 1. The factor (fj,/e)~ 2 in the second term on the right-hand side of Theorem \6.11\ is likely 
suboptimal. 



6.4 Boundary mismatch of the expansion 



Theorem 6.12. There exist constants C , b, 7 > such that under the assumptions 

/i(Af + 1)7 < 1, 

one has 



|W M + Uf L + Vf L || ioo(a/) < C (£) 1/2 e" 



6/m 



Proof. Moreover, at the endpoints of the interval /, the remainder is small. To see this, consider 
the left endpoint of / (the right endpoint is similar). By construction, 



R M (0)|| = U(0)-(W M (0) + Uf L (0) + V^(l/^)) = V A '(1/^) 



tM 



tMi 



(6.14) 
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Theorem 16.81 informs us that Vj\,/ is an entire function. In fact, from Theorem 16.81 and the Sobolev 
embedding theorem in the form IMI^oc^/) — ^II v II.l2(j)IMI.h'1(j) applied to the interval I = [1/fi — 
1,1/ fj] of unit length, allows us to infer 



IvfiCV/i)!!^^ < c||v^|| i2(/) 



I V RlIIl2(/) + || ( v bl) 



< Ce-^\\V% L \\o 



I v rlIIo,/3 + II (Vbl) Ik 



where the constant C > depends only on the choice of /3 made in Theorem 16.81 Recalling the 
definition of v concludes the argument. □ 



6.5 Proof of Theorem 16.11 

The proof of Theorem 16.11 now follows by combining Theorems 16.91 16.111 16.121 with M = 0(l/ji) 
and using the stability result ()2.9|) . 



Acknowledgements: The first author cordially thanks his colleagues W. Auzinger (Vienna) and 
P. Szmolyan (Vienna) for fruitful discussions on the topics the paper. 



A Some miscellaneous results 

Lemma A.l. Let 7 > and 7] S (0, 1). Then the function 

x i-)- f(x) = 7 x (t]x) x , 
is convex on (0, 00) and monotone decreasing on (0,l/(j]je)). 

Proof. We only check the monotonicity assertion. To that end, we compute 

— In fix) = In 7 + 1 + InOtfx) 
dx 

and see that £ lnf(x) < for x G (0, l/{rpfe)). This proves the claim. □ 
The following lemma provides a proof for (|2.1ip and ()2.10p . 

Lemma A.2. Let a > and let BgM 2 be such that £ ■ B|* > a|||j| 2 for all £ G R 2 . Then 

B fcfc > a, k = 1,2, (A.l) 
detB > amax{Bii,B 22 }- (A.2) 

Proof. Property (jA.ip follows immediately from the choice £ = e^, where e/% is the fc-th unit vector. 
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To see property (|A.2p . we start by noting that B 1 is also positive definite: 



x^^x x = y y T B T B- 1 By = y T B T y = (y T B T y) T = y T By > 0. 

Property (|A.2|) follows from properties of the representation of B _1 in terms of the cofactor matrix. 
From 



B 



-i 



1 



detB 



C, C := 



B22 — B12 
— B21 Bn 



and (|A.1|) . we see that detB is positive. The well-known fact that 

MB 1 1| 2 < a' 1 

allows us to conclude 



1 



detB 



|B _1 1| 2 < a' 1 



which implies 



det(B) 



< 



d CI 



Estimating ||C||2 > max{|Cn|, | C22 1 j" and recalling detB > concludes the argument. 



□ 



B Proofs for Section [5] 



Lemma B.l. If u satisfies 



,(») 



00,/ 



<CK n max{n,iT L } n VnG 



for some positive constants C,K independent of (jl, then its complex extension (denoted by u(z)) 
satisfies 

\u(z)\ < CeP^Wi 1 , 
provided dist(z, I) is sufficiently small. 



Proof. Fix x E I and let B r (x) be the ball of radius r centered at x. Then, by Taylor's theorem, 
we have for z £ B r {x) 



\<z)\ < jr 



u^Xx) 



k=0 

LV^J 



A:! 



{z-x) h 



LVa*J 



< 



—k 00 1 k 

y c v_ K n r k + y c * K k r k 

k=0 k=\\/n\ 



< 



1 /M frK\ k 1 00 

l — n \ ™ / , — 1 1 / ,, 1 



k=0 v r ' k=[l/n\ 

If r < l/(2eK) then the second sum above is bounded and we get 

\u{z)\ < Ce rK /» + C < Ce rK/ ». 



□ 



34 



B.l Proof of Lemma 15.21 



We work out the details here: We reduce to the case of homogeneous boundary conditions by 
introducing the linear function u$ with Uq(0) = and uq(1) = <?+. Then the difference u := u — uq 
solves 

~Ijl 2 u" + cu = g — cuq =: g u(0) = u(l) = 

with ||5||z / oo(/) < C g + \g~\ + |<?+|- By the comparison principle (see, e.g., [21 Lemma 2.1] for the 
present form), we have 

„ I. 1 I. I. , , 

< -||<?II-L°°(I)- (B.l) 

From the differential equation, we furthermore get 

// 2 II <> L°° (Z) 2 

\\u < M~ \\g\\L<x>(l) H llsllz,°°(J) < ||p|U°°(/), (B.2) 

where C is a constant that depends solely on c. By an interpolation inequality in Holder spaces 
(see [31 Thm. 3.2.1]) we get from (jB.lj) . (IB.2j) . the estimate 

\\u {1) \\ L oo {I) < c^~ l \\g\\ Lao ^^ (b.3) 

for a constant C > that depends only on c. Combining the above estimates together with 
u = u + uq and recalling that uq is linear, we have 

||« (n) IUoc (/) < Cmax^/i- 1 }", nG {0,1,2}, (B.4) 

where C u = CC g for a constant C that depends solely on the function c. Higher order estimates for 
u are obtained from the differential equation in the standard way. Differentiating the differential 
equation satisfied by u yields 

2„>+2) _ (n) _ A» \ >-!/),.(!/) 



li ' = g 



z/=0 



£ " U-" )«M. (B.5) 



The analyticity of c implies the existence of C c , j c > such that 

||c (n) || L oo (/) < C c7 >! Vn G No- (B.6) 

We claim 

||-u (n) \\ L oo {I) < C u j n max{n, /z" 1 }™ Vn G N . (B.7) 

This estimate is valid for n G {0, 1, 2} by (lB.4j) . To see that it is valid in general, one proceeds by 
induction on n. From (]BT5]) we get in view ()B.6D 



M 2 ||n( n+2 ) 



|£- (J) < ^ 7 n max{n, ii -H n + X;( ri )^7cH!C' c 7 n " l/ max{z,,/i- 1 } n ^ 

i/=0 

n 

< Cgj n max{n, n" 1 }™ + C C C U ^ l v c l n ~ v max{z^, /i" 1 }" 

< <^ 7 n max{n, n" 1 }™ + C C C U max{n, / u" 1 } n 7 n ^ 



y=0 

< C g j n max{n, fi- 1 } 17 + C C C U max{n, ^ 1 } n 7 ™- ' 

< Cu^maxK/z- 1 }"^- 2 ^ 5 f ' 



7c/7 



7 * 7^ + 
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Since the expression in square brackets is bounded by 1 for 7 > 70, if we make 70 sufficiently large, 
we have the desired estimate (|B.7|) for u. 



B.2 Proof of Lemma 15.41 

Set 5 := 25 < 1/e. Define F : k H- 5 k (^Tr,:)* * and note that 

ln(F) = khxS + (i — k) Qn(i + 2) - ln(i + 1 - k)) , 

i — k 



(ln(F))' = hx5-hx{i + 2) + hx{i + l-k) + 



i + l-k 



< hit + In ( 2 k ) + 1 1 I + hi'l 



Since i < 1/e we have (ln(-F)) < 0, and F is monotone decreasing which gives 

fc=0 V 7 A:=0 



B.3 Proof of Theorem 1531 



We will only consider the inner expansions , at the left endpoint of /, since for the expansions 
at the right endpoint the arguments are almost identical. The proof is by induction on i. For i = 
we have from (|5.4j) - (|5.7j) . 

f _/,2„« 1 (^22011^012021) _ _ 021 f 

) V V + an U ° ~ 9 an / (B 8) 

\ «o(0)=«„(l)=0 ' 1 ' ^ 

u = — (/ - ai2^o) , (B.9) 
an 

-(«§-)" + on (0)t^ = 
u£(0) = -«o(0) 



(B.10) 

/"OO /*OG 

v% = - / a 2 i(0)u^(T)dTdt. (B.ll) 



Combining Lemmas 15.31 and 15. 2| we can find a constant C4 > such that satisfies (j5. 13j) . which 
in turn shows that uq satisfies (|5.12p by Lemma 15.31 The solution formula for Uq then shows that 
(|5.1ip is valid. Finally, Lemma [4T61 establishes the desired bound (|5.10p for v 2 - 

So, assume (|5.10p ^ (|5.13|) hold for up to i > and establish them for i + 1. We will choose the 
constants Oi such that the ratios C2/C1, C3/C2, and C4/C3 are sufficiently small. Furthermore, 
the constants Ki and C u , C v are sufficiently large and are selected such that 

Ki=K 2 K 3 = K 4 , C U = C V (B.12a) 

C%K 1 = K 3 C u = a>2/min{l/a 11 (0) ! l/ai 1 (l)}. (B.12b) 
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We start by bounding and Ui + \. We first consider u^+i, which satisfies (see eq. (|5.4p ): 



( „" _i_ ("22011-012^21) -, _ ,,2 021 ,,// 

\ «,-+i(0) = -vf +1 (0) , v i+1 (l) = -vf +1 (0) 

To that end, we assume, as we may, that 7 is sufficiently large for Lemma 15.31 to be applicable 
for the function Ui-%. Together with the induction hypothesis, it then yields Vn E No (with the 
constant C of Lemma |5, 31 for the function g = 021/011) 



2 ( °2i // 

1 V 

From 



(n) 



If — u'U) || L oo (/) < ^C'C 3 K l 3 -'((i + l + n)n + 



(i-l)l 



7 n+2 max{n + 2,^ i } 



li n+2 



n + 2 



n+2 



max{(n + 2)//, l}* 



/j 2 max{n + 2,/i" 1 }" +2 < max{n, /i _1 } n 

< C"max{n,/x- 1 } n ((i + n+l)/x + l) 2 
for some C" > independent of n, we get Vn E No 

( n ) f„- _ 1 v-i 



II/* 2 H*"M ^ CrfC'C"K^K^((i + l + n)(,+ 

The induction hypothesis for vh_ 1 and and Lemma [5.2l therefore produce (with C of Lemma l5.2 
for all n E N 



(1 - 1)! 



lin 



i II L»m < C7 n max{n, /i } 



fi-1 



(i-l)I 



i-1 



CbiQ-^Ci + 1 + n)n + l) i+ V + 2dK\-' + — - ) (i - l)^ 1 )^*" 1 ) 
< C 4 Kl +1 -f n max{n, /i _1 } n (/i(i + n) + if 



(j _|_ 1 
(i + 1)! 

C^cg7 2 ^" 2 (^3/K4) i+1 + KMC ' g(CX/^4) m 



< C 4 i^ + V max{n, /i~T(M(* + n) + 1)* 



^ ) %^ 2 + 2C(-)c(^) : 



where we made use of - rez^ri < C (1+IT! anc ^ ^ e si m plify m g assumptions on the relations 
between the constants K{ in (|B.12|h For the induction argument to work, we have to require that 
the expression in square brackets is bounded by 1: 



< 1. 



(B.13) 



We will collect further conditions on the constants and Ki and see at the end of the proof that 
these conditions can be met. 



37 



Next, we turn to to u i+ i, which is given by ([5.50 : 

1 



Ui+l(z) 



an (2) 



[/u 2 <_i(z) + 012(2)^+1(2;)] 



Lemma 15731 the induction hypothesis and the just proved bounds for i> i+ i then give for a constant 
C that depends solely on the coefficients of the differential operator, 



l4+ilU~(l) < C 3 ^ +1 ( / u(z + l + n) + l) i+1 max{n,^- 1 r ( ' + l) " 1 



(i + iy. 



In view of K3 = K4 by (|B.12p . we recognize a second condition for the induction argument, namely, 



^3 



< 1. 



(B.14) 



Next, we consider uk_ l which satisfies (see (j5.6|) ): 



{uf +1 )" + an(0)uf +1 = -a 12 (0)Sf +1 - £ 



i+l / (fc) /„n 

k \ A 4 "j+l-fc "T fc] P X "i+l-fc 



fc=l 



a i2 } ( ) ,,k~k~L 



uf +1 (0) = -u i+ i(0) , -+ as x -+ 00 

(B.15) 

We will estimate the right-hand side of (|B.15p . To that end, since 1 < k < i + 1, we start with the 
observation 



1" A* + 



1 



i + 2-k 



i+l— A 



< A*" M + 



1 \ 



i+l— fe 



i + 2 



< A* 



i + 2 



i+l 



i + 2 



i + 2 \ 
j + 2 - aJ 

i+l-fc 



i+l— k 



i + 2-k 



so that we can estimate 



it=i 



1 v +1 i 



i + 2 



(C^ + l) + |z| 



\2i+lg-/9Re(z) 



i+l / \ fe 

2±] 



i + 2 



i+l-k 



If is sufficiently large, then the sum is bounded by 2-^e by Lemma [5741 Analogously, we get 
(note that the case i = leads to empty sums) 



k\ |z| M |Ui+1 - fcl ^ k\ 

k=0 k=0 



gl«S?(0), 



PI A 4 |«»+i-aI 



< c.Ci^r 1 („ + i)" 1 - 1) + Mf-V^) £ 



fc=0 



7a 



k / \ i—l—k 

I 



i — k 



<K^C"'C a C l K{^[^ + 



i+l 



n i+l 1 



(C v (i + l) + \z 



\2i+l -0Re(z) 
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where we appealed again to Lemma 15.41 and used that K\ = K 2 is sufficiently large to bound the 
sum by 2e, and noticed additionally that 



1 



i-1 



< U t + 



1 



i + 2 



i-i + r i-i 



< M + 



1 



i + 2 



t+i 



2 + 1) SUp 



i>l 



i + 1 



=:C""/(2e) 



where the last supremum is finite. 

These estimates allow us to bound uf +l with the aid of Lemma 14.51 to arrive at 



\ut +1 (z)\ < C 2 K* 1 [ l i + 



1 



i+l 



-/3Re(z)_ 



1 



i + 2 
2 Ci /i^i 



(C v (i + l) + \z\ 



\2(i+l) 



(i + l)! 

1 c 2 \k 2 ) + c 2 \c 2 v k 1 



Together with the simplifying assumptions (|B.12p we see that for the induction argument to work, 
we need to require 

C^ + K-^ + ^]<l. (B.16) 

K 2 O2 O2 



Finally, for u£_ 3 we have from (|5.7p 



^+ 3 (*)l 



E[ Q 21 (°) ,,fc^fc-L , a 22 ; (°) 
M ^ T U i+l-fc + 



fe=0 



A:! 



dfdt 



< ^Imaxl^CO)!, |4?(0)|}J / f? k uf_ k d?dt + / [r k vt k d?dt 

k=0 [ z t z t 



Proceeding analogously as before, we find that this last sum can be bounded as 



\rhs{z)\ < C 1 K[ +1 [n + 



1 



i + 2 



i+l 



(i + l)! 



(C v (i + l) + \z\ 



^2(i+l) -/3Re(z) 



c 



C 2 (K 2 



i+l 



+ Ci^f 



where the constant C > is suitably chosen. Lemma 14.61 then gives 



\vU*)\ < ^ +1 (^ + ^) 4+1 ^^(a(. + l) + k|) 2 ( i+1 )e-^^) 



c 9i (Zi 

Ci V-Ki 



i+l 



/? 2 Vl-(2(i + l))/(/9C„(t + l)) 



Hence, by our assumption (jB.12|) . we see that we can find C (depending solely on 0) such that our 
induction argument will work if we can satisfy 



a 



C- 2 - + CK{ 2 
C>1 



< 1. 



(B.17) 



In total, we have completed the induction argument if we can select the constants Cj and Ki such 
that (|B.13|) . (|B.14p . (|B.16p . and (|B.17p are satisfied. Inspection shows that this is the case by 
taking the K sufficiently large and appopriately controlling the ratios C 2 /C\, C3/C2, and C4/C3. 



39 



B.4 Proof of Theorem 15771 



From Lemma 15.61 and the estimates (|4.2ip , (|4.22h , we get 

M 



fjAf 



(n) 



x) 



< Tr > i - ! !//, -h Ji'-; ••' C< > I Me + — I 7i < 

i=0 i=0 ' 

where, in the last step we have used the assumption that e(M + 1 + e/fJ>) is sufficiently small so 
that the sum can be estimated by a convergent geometric series. 



B.5 Proof of Lemma 15.101 

We start by noting that the change of summation index k to t = i — k leads us to having to estimate 

i>M+l VP/ £=0 V 7 

We start with the elementary observation (cf. Lemma I4.8P 



for suitable 7 > 0. Next, we estimate 

where we suitably chose 7 independent of £ and 5?. We therefore conclude 

M 



E ^)'E" , -'^r'(M+^ I ) i(w+s)*.H- 

1_A V /v 



£=0 

„ M 
< _^_ e -/3£/2 £ M+l + ^ 

1 — ^— 

£=0 

Again, elementary considerations (cf. Lemma HTH]) show x M+1 ~ e e~ l3x ^ 4: < r )^ I+l ~ l {M + 1 — £) M+1 ~ 
for suitable 71 > 0, so that we arrive at 

M 

S < -^ e -^/2 e ^+iV(7A7l) M+1 ^(M + l-^) M+1 -Y(^+l/^ 

1 — — J 

e=o 

< _E^_ e -^/2 £ M+l 7 M+l (M + j + 

1 — X 

where, in the last step we have used the convexity of the function 

i a M+1 -"{M + 1 - £) M+1 - £ (£ + 1/ M )< 
on the intervall [0, M] - see also the related Lemma ID. II 
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C Proofs for Section [6] 



C.l Proof of Lemma 16.41 



The proof consists of showing that the terms of the sum can be bounded by the terms of a convergent 
geometric series if 5 is sufficiently small. To simplify some notation, we restrict our attention to 
the case i > 2 (the cases i = and i = 1 being easily seen). We denote the terms of the sum by 
F(n) and note that it suffices to study the sum Y^n=i F{n) since the case n = produces the term 



(r + m) 4 (rre + l) m+1 \ r + rre / \ m + 1 

which is bounded uniformly in i if 5 < 5q < e~ l . 

We define, for re G [l,i — 1], 
/(n) := lnF(re) 



and compute 



(r — 1 — n) In (5 + n ln(n + m) + (i + n + 1 + m) ln(r + re + 1 + rre) 
— r ln(r + m) — (2re + 1 + rre) ln(2n + 1 + m) 



, (n + m)(n + i + 1 + m) n 
/ (re) = In In d — 1 + 



(2re + 1 + m) 2 re + m 

It is easy to see that f'(n) > —1 — In 4 — In 5 for n > 1 so that we can find a constant c > with 
f'(n) > c > for all re > 1 by taking 5 sufficiently small. By the mean value theorem, we therefore 
get for each re, that In — fejy^ = /(re + 1) — /(re) > c. Hence, i^n) < e~ c F(n + 1). Iterating this 
estimate, we get 

F(re) < e- {i - l ~ n)c F(i - 1). 

Hence, 



i-1 i-1 



E F ^ ^ E e-^ 1 -^^ - 1) < 73^* " !)• 

n=l n=l 

The argument is concluded by noting that F(i— 1) is bounded uniformly in « and rre as the following 
rearrangement shows: 

(i - 1 + rre) i-1 (2i + rei) 2i+m _ 2i + m /i - 1 + m V' 1 ( 2i + m \ 2i + m ~ 1 
^ '~ {i + mf (2i + rre - l)2*+m-i ~~ £ + m ^ i + m / \2z + rri — 1 

C . 2 Proof of Theorem IBTTl 

Specifically, differentiating rre times yields (we write Uj instead of Uf) with B = A(0) 

i-1 



E u(-+2) = Bu; m )+E(y- n A l _ n u n ) (m) 

n=0 

i-1 min{m,i— ra} / \ / ■ \ 

= bu^ + Ea_ E ("•')( / "V' r ' " ' 

n=0 j=0 V J / V J / 
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We abbreviate 

min{m,i— n} . . / . 

F n := g (7)( , 7)iS i - W -^" ) , 
and we proceed as in the proof of Theorem l6.61 The induction hypothesis yields, for any f3 S /3 t 



oj 



x>0 



(• \ i—n—j 
1 - U ~ J (fa - /3)-( 2n+1+m )(n + m - j) n K?K!T j v-<> m -fi. 

(3 ~ P J 

For j <i-n,we select /3 = f3 + re(/3 - P) with k = and get 

\\$-"-i\j<F-i)\\ Qtf} < (C.2) 
• , n (i 4- n — i 4- 1 4- m ) i + n ~3+ 1 + m 
"V" (2n + 1 + m) 2n+1+m V J7 1 2 

the induction hypothesis shows that the estimate (|C.2|) is also true for j = i — n. Therefore, by 
estimating n + m — j < n 4- m, we get 

||F n || 0j/3 < Ca^o-/?)-^ 1 ^^?^!/^- (/, + '" ) " 



(2n + 1 + m) 2 ^^" 1 

min{i— n,m} 

X 



E (7) C i y i!(A} " /i,A ' 2 ' v,{; + n " j + 1 + ™r +n ~ j+1+m 

With the estimates ( m )j! < m J we bound 

min{i— n,m} 



I — "5 " J / \ /" \ 

e (7) C j n ) _ ^ Vi/ + n - 3 + 1 + ™) t+rt ^ +i+m 



\i+n+l+m 



3=0 



= (1 + (A) - P)vK^ l y "(i + n+l+m) 
where we have recognized a binomial sum in the last equality. Hence, 



i-l 



^HA^yFjo,/, < C u C AlA {fa-P)-^ +1+ ^K{- 1 K?v- 



n=0 

i-l 



X 



^ x _ m , iy - n (n + mni + n + 1 + m )i+n+l+m 



n=0 



(2n + 1 + m )2n+l+m 

< C^CoCataGSo - /?)- {2l+1+m) Kr 1 K 2 m ^- m (i + m)\ 

where we appealed to Lemma 16,41 and implicitly used that K\ and Ki are such that ^ A K^(\ 4 
(/3o — ^vK^ 1 ) is sufficiently small. Using 1 1 IE1 1 1 1 2 < v~ 2 , we therefore get from (jC.ip 

||ul m+2) || 0l/3 < ^ 2 Q/(A) - /3)-^ +1+m )(i + m)^l^ 2 m ^" m [l|B|| 2 + C C A1A K^] 

< Cu((3 - (3)-^ +1+m \i + m) i K{K^ +2 u^ m+ ^ [i^ 2 ||B|| 2 + K^ l C C AlA K^ 1 } ; 

the expression in square brackets can be bounded by 1 if K\ and Ki are sufficiently large. 
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D Proofs for Section |4] 



D.l Proof of Lemma 14.21 



The proof is by induction on i. For i = 0, equation (|4.17p gives all the assertions (trivially). So, 
assume (|4.18p - (|4.2Up hold and establish them for i + 1. 



We first consider (|4.18p and assuming 



H-2,j-2 
V i-2,j 



we want to show that 



u i+l,j 
Vi+l,j 



t-i,i-2 



V j >i, 



V j > i + 1. 



By the induction hypothesis and j > ? + 1 > i — 1, we have Vi-i,j = 0. Also, from j > i + 1 we 
get j — 2 = j — 1 — 1 > i — 1, so that the induction hypothesis implies Uj_ij_2 = 0. Therefore, the 
right-hand side of (|4.17j) vanishes and thus the induction step for (|4.18p is accomplished. 



For (|4.19p let :={(i,j) £ZxZ: t<I, i odd or j odd}. We proceed by induction on I by 

assuming (I4.19|) to be true up toX. For i < I with G we have that G A4(T+1) 

implies (i — 1, j) G M(I) (either i + 1 is odd and then i — 1 is odd or j is odd) and additionally 
(i — 1, j — 2) G Hence, for (i + 1, j) 6M(1+ 1), by the induction hypothesis, the right-hand 

side of (|4.17p vanishes, which proves the induction step. 



We finally consider (|4.20p and we want to show 

\u i+hj (z)\ + \v i+1J (z)\ < C s r i - 1 K i+1 {i + VzeGj. 

We set 

C A ■■= sup||A _1 (x)|| £ i, 

where ||(xi, X2)||^i := |xi| + |x2| denotes the usual t -norm and ensure that K satisfies 2C A /K 2 < 1. 



Let k G (0,1). By (14. 17ft . the induction hypothesis with G^^^s C Gs and Cauchy's Integral 
Theorem, we have 

+ k+ijO)! < Ca (K_ij_ 2 (z)| + K-ij(«)|) 

1 1 2C A 



(k5) 2 

< C s S~ i - 1 K i+1 (i + l) ( - i+1 ^> 



K 2 (i + 1) 2 k 2 (1- nf- 1 + l 



i-1 



The choice k = gives 



\u i+ltj (z)\ + \v i+1>j (z)\ < CVT^iT+^ + l)^ 1 ) 
by the choice of K. 



C A 
2K 2 



<Cs5~ {i+1) K l+1 (i + l) {i+1 *> 
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D.2 Proof of Lemma 14.51 



We provide some details for the case Rez < 0. For z G (0, oo), the use of a Green's function gives 
the following representation of the solution u(z): 



u(z) 



1 

2^2 



e y f(y/a)dy + e az e-*f(y/a)dy 

J a 



az 
az 



Analytic continuation then removes the restriction to (0, oo). In order to get the desired bound, we 
estimate each of these four terms separately. We restrict here our attention to the case of Rez < 0. 

For the first integral, we use as the path of integration the straight line connecting and az to get 
e y f{y/a)dy 



< e - Re{az) f C f (q + t\z\) j \az\e- Re{taz) e Re{taz) dt 
Jo 

< e - Re(Tiz) [ C f (q + t\z\y\az\e Re{{1 - t)z(7I - a)) dt 

Jo 

< c f e- R ^-^- i {( q + \z\y +1 - q >+ 1 }. 



For the third integral, we calculate with [H eq. 8.353.5] and the incomplete Gamma-function •), 



e y f(y/a)dy 



< C f 



-y-a/ay 



{q + y/a) 3 dy 



= C f a- j {l + a/«) _(j+1) e a9(1+ - /a) r(j + 1, aq(l + a/a)) 
= C f a(a + a)- (i+ V (a+ ^r(j + 1, q(a + a)). 

In view of the assumption (a + a)q > 2j + 1 > j, we may employ the estimate 



|r(a,OI < 



iei-ao' 

(see, e.g., [121 Chap. 4, Sec. 10]) to arrive at 



«o = max {a — 1,0}, Re(£) > 0, |^| > ao, 



e~ y f(y/a) 

Hence, the third integral can be estimated by 



< C f aqi +1 



1 



q{a + a)-j 



(D.l) 



e- y f(y/a) 



< C f aq 1+l 



1 



-a Re(2) 



q(a + a)-j 

We now turn to the second integral in the representation formula for u. We split the integral as 



e y f(y/a)dy 





raz 


< 


az 1 

e e 




Jo 



e v f{y/a)dy 



+ 


e az [ 




Jo 



e y f(y/a)dy 
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We recognize that the second integral can be estimated using (jD.ip . The first integral is treated as 
follows: 



''■"'/ e y f(y/a)dy 
o 



< C f e Re{az) I e- tRe{az) \az\e- tRc( - za \q + t\z\y dt 



< C f e- aRe{z) f 1 {q + t\z\) j a\z\e^- t ^ a+a ^ Rez dt 

Jo 

< C f e- uRe ^— Uq + \z\y +1 - q> +1 ] . 

j + 1 L J 

Combining the above estimates and recalling q(a + a) — j > 2j + 1 — j > j + 1, we conclude 



raz poo roc 

' az / e y f(y/a) dy + e az / e^f(y/a) dy - e~ az / e -»/(y/a) dy 

JO Jaz JO 

a 



< 



2C f e- Re(az) (q + \z\) j+1 -^ 



3 + 

Combining this estimate with the obvious one for the fourth term, we arrive at the desired bound. 



D.3 Proof of Theorem WM 



We begin by setting 



a r- 



det A(0) aii(0)a 22 (0) - ai 2 (0)a 2 i(0) 



au(0) a u (0) 
and choose the constants C„, CV, C„, C^, Cj, Ki,Ki, i = 1, 4 to satisfy the following: 

K 1 = K 2 = K 3 = K 4 > 1, 
TTi = K 2 = if 3 = i? 4 > l, 
C« = = C„ = > 1. 

Furthermore, the following requirements have to be satisfied: with the constants Cs, K of Lemma l4.2 
(we assume for notational simplicity that 5 = 1 is admissible in Lemma l4.2p and 7 a the constant 
of analyticity of the data (see (|2.3p ): 



(D.2a) 
(D.2b) 
(D.2c) 
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C 2 > Cs, 
C 3 > C S , 

Ol2 



Ci > 
Ci > 

c 4 > 



c 2 , 
c 2 , 

2 C 3, 



on(0) 

|02l(Q)| 

a 22 (0) 
|o 2 i(0)| 



Ki = K 2 > max{l,7 }, 
1 > 



1 > 



1 > 



1 > 



1 > 



> 



1 > 



1 > 
1 > 



cJ 1 + ^MU( 1 + £ 

V an(0) J a' V C 2 

g 2 |ai2(0)| C a 7 a 1 f 1 C 2 
d an(0) + a n (0) K x 1 - 7 a/#i V Ci 



1 la C s (K 
l- la /K 2 K 2 C 2 \K 2 



i+l 



c s ( i y +1 , a 



C3 V^3 



C 3 |o 21 (0)| 



1? \ J+ 1 



C 3 V^3 



+ Kf^\a 12 (0)\ 

^3 



C 4 a 2 1 - 2/(C e a) 



Cs ftf \ 



i+l 



c 3 V^s7 



1 



1 - 2/(Cea) 



C a 



C3 + 7T 2 

c^ 4 



i 



1 - la /(K 4 K 4 ) 



C 3 , 4—2 



1 > ^1021(0)1^ + ^^4 022(0), 



C 4 a 2 

2Ci 2e»| Q2 i(0)| Cs_ H V +1 r 2^ /^EiV +1 
C 2 a'a 22 (0) C 2 \kJ 4 C 2 V^J 



|ai 2 (0)| C 2 __2 2e a 
+ 



an(0) Ci 



on(0) 



7a 



021(0)1 g Cx z -2 / [021(0)1 

L «n(0) le C 2 2 + ° a V a n (0) ; tfi 1 - 7 a/*i VC 2 



c 2 v^j 



lv+ wg(£ 



£4y +1 

K 2 J 



Oi2(0)| C 2 _- 2 2e« 

+ 1 ~ + 



1 



7a 



on(0) Ci 
1 



,1 + ^ 

a n (0) an(0)l- 7a /^iK 1 V Ci 



on(0) 



7a 



C a 



C 9±R- 2 + 
a C 3 4 K 3 ^ 3 1 - la /(K 3 K 3 ) 



+ 



Cs f K 



C 3 V^3 



1 



#3 



+ 



Ci (K 



-i 



C 3 V^3/ V^S 



^3 



(D.3) 
(D.4) 

(D.5) 

(D.6) 

(D.7) 
(D.8) 

(D.9) 
(D.10) 

(D.ll) 

(D.12) 
(D.13) 

(D.14) 

(D.15) 
(D.16) 

(D.17) 
(D.18) 
(D.19) 

(D.20) 
(D.21) 
(D.22) 
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Before we proceed with the proof, we make sure that the requirements (|D.3|) ()D.22|) can be satisfied. 
We make the following simplifying assumptions in addition to (|D.2P : 



C u 



C v 



C 2 = Q, 



la _ la 

K x ~ K 2 
K <K\ 

C 3 = 



— C u - 

_ la_ 

K 3 
= K 2 



Gs > 8/ a, 
= 2^- 

: K 3 = K 4 , 
Ci = Q A . 



1 



(D.23a) 

(D.23b) 

(D.23c) 
(D.23d) 



Here, Q > will be selected sufficiently large below. Then, the requirements (jD.3h (jD.22h are 
satisfied if: 



Q > C s , 



Q 3 
Q 2 

Q 2 

Q 4 



> C s , 

l«12(0)| 

aii(O) 

l«2l(0)| 



> 



> 



> 



1 > 



1 > 



1 > 



1 > 



1 > 



1 > 



1 > 



1 > 



1 > 



1 > 



1 > 



a 22 (0) 
I "21(0)1 



Q, 
Q, 



C a 1 + 



|Q2l(0)| 

au(0) 



2 
a 7 



Q |ai 2 (0)| C a 2i a 
Q 2 ou(0) au(0) ifi 

Q 3 2|a 2 i(0)| 
Q 4 a 2 

g 3 ' 

4 



Q 2 


\ 7a 


+^1 

Q. 


Q 


J K 2 


i + 






Q 2 J 


5 



2C„ 



4 



g 4+J ^ 4 . 



,<9 , 1—2 

a 2 



1 > ^21(0)1^4+^4 «22(0), 



_- 2 g 2 2e> 2 i_ 
1 g a'a 22 (0) 



+ Q +Ka Q' 



|ai 2 (0)| g ^-2 2e Q 
au(O) g 2 1 an(0) 

oii(o) g 



ai 2 (0)| g , -^-2 2e a 



au(o) g 2 
1 



an(0) 

2Cg 7 a 

au(0) ' an(0)Ki 



021(0)1^ 2 7a fg 2 

g 



A'i 



i + i ( — - 1 1 + — a r^- 



g 



+ 



i + 



Q_ 

Q 2 



au(0) 



J7a_ 



g 3 g 3 
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It is now easy to see that by first selecting Q sufficiently large and then choosing the parameters 
Ki and Ki sufficiently large, ensures the above requirements. 



We now turn to the induction argument. 



(Base) Case i = j = 0: 

The functions Uq ,v^ ,Uq ,Vq 2 C?o,o = 0) satisfy the following: 

J _ -Ro)" + «% = 

\ v o,o(°) = -«o,o(0) , ^ (x) -)• as x ->■ oo 

-L «12(0)~i 



(D.41) 



"o,o = -TTW^o, (D.42) 
an(U) 



J -( S o,o) + a n(0)^ = (D43 x 

\ S o,o(°) = - u o,o(0) , «o,o(^) ^Oasx^oo ' 

roo poo 

Vo, 2 (z)= / / a 21 (0)u^ (r)dTdt. (D.44) 

./z 

Solution formulas for Vq and Mq and Lemma I3~2l (recall that we assume that 5 = 1 is admissible 
in Lemma [42]) give us the desired result for Vq and Uq in view of the requirements (|D.3|) . (|D.4p . 
while (|D.42p gives it for Uq , in view of requirement ()D.6|) . For Vq 2 (z) we have from (jD.44|) . the 
just proven result for Uq , and Lemma [3~6l (with j = 0) 

\vL(z)\ < |a 2 i(0)|C 3 ^Exp(z) < C 4 Exp(z), 

by the choice of C4 in ()D.7p . 
(Base) Case j = 0, i > 0: 

The functions uf ,vf ,uf ,v^ 2 satisfy equations (|4.23p - (|4.27p . respectively, for up to i. We proceed 
by induction on i. First, vh 1 satisfies (14.24p with i replaced by i + 1. The right hand side of that 
boundary value problem satisfies, in view of the induction hypothesis and the choices K\ = K 2 , 
Cu = Cy > 1, 

^( 1 + ^)W(«TTr^( c W i + 1 -*) + l' 



2(i+l-fc) x 



fe=l 

x (CiA'j +1 - fc + C 2 ^ +1_fe ) Exp(z) 



l Q 2l(0)l \ r n _ , \ 1 7a 1 . , I i\2t+l 



< C a 1 + (d + C 2 )- JS-Kf^Qfl + Exp(z), 

aii(0) y 1 - 7a/K 2 A 2 i! 
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where we used the fact that K\ = K<i > j a by (|D.8|) . Lemma 14.51 and Lemma 14.21 yield, for the 
solution v-'q of the boundary value problem (|4.24|) . 

+ C s K l+ \i + iy +1 )Exp(z) 



c a 1 + 



t»2i(0)h i A , Ci 



i + 



I 7a C 5 (K 



a u (0) ; a' V l-7a/^2^2 C 2 



i+l 



(D.45) 



The expression in square brackets is bounded by 1 by our requirement (|D.9p . Next, we consider 
Ui+i o which satisfies (|4.23p with i replaced by % + 1. We have by the induction hypothesis and 



\u. 



+1,0 



Ml < { 



|«12 (0)| 



(TPi)! (CV(i + 1} + kl)2(i+1) + 

Ci . fc I ife 1 



au(0) 

i+l 



fc -a n (0) ' al ~' (i + l-fc)! 
+ C 2 ^ +1 " fc (C 7 (i + 1) + |z|) 2 ( i+1 - fc ) 



Ciiq +1 -* (<%(« + 1) + \z\) 2{l+1 - k) 

} Exp(z) 



C 2 ki2 



+ 



C a 7a 1 



d au(0) au(0) ifi 1 - 7a /ifi V C x 



1 + 



Co 



Exp(z) 



since K\ = i<C 2 > 7a by (ID.8j) . Again, the expression in square brackets is bounded by 1 by our 
requirement (|D,10p , 



For uf +l , we have by Lemma 1431 (|4.26p . and Lemma I4~2 



\uf+i,o(z)\ < e 



-an(O)ReOz) 



-i,o(0)| <C s K l+l (i + lf + V Exp(z) 



C s ( K 



i+l 



Exp(z) 



(D.46) 



c 3 V^3 

In view of our requirement (|D.14p . the expression in square brackets is bounded by 1. 

Finally, for v^ +l 2 we have from (|4.27p . ()D.46P and Lemma I4~6l in view of K 3 = and = Q. 



1^+1,2^)1 < 



a2i(0)uf +10 (r)dTdt 



z t 



< 



|<a 2 i (0)| 



a 2 \l-2/(C%a 
< C 4 Kl +1 (Cv(i + 1) + |z|) 2 ( i+1 ) Exp(z) 



C 3 Kl +1 (Cz(i + l) + \z\) 2 ^ Exp(; 
C 3 |a 2 i(0)| 1 



C 4 a 2 1 - 2/(C 2 a) 
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The expression in square brackets is bounded by 1 by our requirements (|D.13p . (|D.12p . 
(Base) Case i = 0, j > 0: 

The functions j+i'^o j+i'^o j+i'^0j+3 satisfy the following: 



o,j+i ) + a ' v o,j+i 



»2i (Q) (~L 



W) 1 Oj-i 



oj+i(°) = " Kj+i(0) + "o,i+i(°)) . v o,j+i( x ) ^ as x ^ oo 



U 0J+l 



' V 0,j+1 



ai 2 (0)~ L l u o,j-i 



an(0) 



ail(0) 

(ufe+l) +an(0)^ J+1 = -a 12 (0)^ J+1 



(0) = - U j+i(0) + «o,i+i(°) ) ' S o,j ^ as x ^ oo 



Wo,i+3(^)= / / [fl2i (0)^+1 (r) +a 22 (0)i5oj+i(r)] drdt. 



(D.47) 

(D.48) 
(D.49) 

(D.50) 



2 t 



To establish the desired claims, we proceed by induction on j noting that the case j = and i = 
has been proved already. Assuming that the bounds are valid for i = and up to j, we show them 
for j + 1. 

We start wit h Vqj + 3 - By the induction hypothesis, Lemma l4~6| and the assumption Cga > 4 (see 
requirement (|D.12|) ) we have 



< ^ (\a 2 i(0)\C 3 Ki + \c q (j + 1) + 2tt+D^ 



0' + 1)1 



+ a 22 (0)C 4 ^i \Cv{j - 1) + |^|)20-i) 1. 



(i-i)! 



Exp(z) 



In view of = and = K4 we get 



f', 



0,j+3 



'z)\ < C A K 



4^4 



-(C^ + l) + |z|)^ +1 ) Exp(z 



4 4 2 

-|«2i(0)|-f + -|a 22 (0)|K 4 

(2 O4 tt 



(i + i)! 

by requirement (|D.16j) . the expression in square brackets is bounded by 1 as required. 
We next turn our attention to «o,i+i wmc h satisfies (jD.47|) and right-hand side (RHS) satisfying 

|O2l(0)| 



RHS HDA7\ 



< 



\an\ 



(D.51) 



We first study the case j = 0. Then Uoj_\ = 0, and Lemma [431 yields, together with Lemma I4T2 

| e -a'Re(*) Cs < c 2 K j 2 +l Exp(z) 



l^oj+ifc) 



c^ 2 
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Since we assume K 2 > 1, our requirement ()D.3p implies the desired bound. Returning to (|D.5ip 
for the case j > 1, the induction hypothesis and Lemma 14.71 produce 



RHS mM (z)\ < lMM 2e ^ 1 ^- y mj - 1) + \z\ + lf {j - l) Exp(z). 
Therefore, Lemma 14.51 together with Lemma 14.21 yields 



< 



73-1 



1 1 o „>2i(0)L K{ (Cfe (,--l) + k| + l)^ +1 + |«,y + i(0)| + |^ +1 (0)| 



-2e 



Ci- 



a'2j-l a 22 (0) 



Exp(z) 



< c^ +1 — + 1) + \z\f^ 



a +1)1 



2 Ci2e> 21 



C, ( 1 



i+i 



C 2 \Kx 



Exp(z), 



C 2 a'a 22 (0) C 2 V^i 
where, again, the expression square brackets is bounded by 1 due to our requirement (|D.17p . 
Now, consider j+i' w hich satisfies (|D.48p . By ()D.52p and Lemma UTT] we have 

l«12(0)| 



(D.52) 



< 



on(0) 



■C2K ] 2 +1 TT ^(C~ v {j + l) + \z\\^ 



(7 + 1)1 



2e a ^^,-1 1 



an(0) 



(i - 1)! 
(C«0" + 1) + I* 



(Cat? " 1) + 1 + M 



\2(j-l) 



Exp(z) 



< g^* 1 — L_iY' 1 > + 1 ^ + l-h 2 "+" 



|ai 2 (0)| C 2 , =-2 2e a 



on(O) Ci 



an(0) 



Exp(z); 



in view of requirement ()D.18p . the expression in square brackets is bounded by 1 as required . 

Finally, for Uq - +1 which satisfies (|D.49P we have by Lemma 14.51 (the case j = needs special 
treatment in that the third term in the following estimate is not present) 



17 j~ l 



ko,i+i(0)| + K J+ i(0)| + \a 12 (0)\C 4 Ki 



< c^j-^mi + 1) + \ z \^) 



1 



(j - l)!(2(j - 1) + 1) 



(C d (j + l) + \z\ 



(i + i)! 

Again, in view of our requirement (jD.lip . the expression in square brackets is bounded by 1. 



Induction Step: We proceed by induction on j, the induction hypothesis being that the estimates 
(|4.32p - (j4.35p have been shown for all i G No up to j and will establish them for all i £ No and j + 



Exp(z) 



Exp(z). 



We start with estimating - +3 . In view of the definition (|4.3ip . we estimate with the induction 
hypothesis (recall K3 = K4 and K% = K4; we also point out that for the case j = 0, the terms 
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stemming from j+i—k are i n fact not present) 



min{i,j+l} fe 



k=0 



fc! 



min{j j'+l} 



< 



^ E ^ 



k=0 



C 3 K l f k K j 3 +1 k 



+ j + 2k) 



-(Cz(i + j + l-2k) + \z\) 2 ( i 



+i+i-2fc) 



(i + j-l- 2k)\ 



(C q (i + j-l-2k) + \z\) 2 ^ 



i+j-2k-l) 



Exp(z) 



<c 4 ^ +1 ^ + ^ + 1 ) + l^ 2(i+i+1) 



c 4 + 4 



1 - Ja/(K A K 4 ) 



here, we employed observations of the form (note that i + j — 2k — 1 > 0) 



=— Exp(z); 



< 



(i + j-l -2fc)! 
(i + j + l) 2fc+2 

(i + j + 1)! 
(i + j + l) 2fc + 2 
(i + j + 1)! 



(Cfc(t + j - 1 - 2k) + \z\) 2( - i+j ~ 2k -V < 



(i + j-l- 2k)\ 



(Cz(i + j + l) + \z\) 2 ^- 



2k- 



(Cz(i+j + l) + \z\) 2 ( i+ i- 2k -V 

(c s (i + j + 1) + \z\r 2k - 2 {c q {i + j + 1) + izD 2 ^'- fc ) 



< + + i + 1) + N) 2(i+J) < (iT7+i)! ( ^ + j + + l^) 2(i+J+1) - 

From (|4.3ip and Lemma 14.61 we therefore get 

|«&+3(*)| < C4KIK4 -, , I ,m Exp(z) 



a 2 Vl-2/(C^a 
By requirement (|D.15p . this is the desired estimate 



i + j + 1)! 

2 

C a 



O4 



1 - j a /(K A K A ) 



We now consider v^j +1 which satisfies (14,29ap -( l4.29bp . The required estimate is proved by induction 
in i, the case i = having been studied previously. The right-hand side of the boundary value 
problem satisfies 



RHS K2M - K2%\ 



< 



l«2l(0)| 



oii(0) 

The induction hypothesis and Lemma 14.71 produce 
RHS {ZM^ - K2lM 



tj-lfW +Ca (l + l -^§) £#1*1* {\ut W (z)\ + \vtk, j+ l( 



< C 2 -^— (C a (i + j) + \z\) 2 ^ +1 K\K{ +1 Exp(z) 

(* + jy- 



au(0) 



2e a K 2 ^ + C a ( 1 + 



a 2 i(0)h 7a 1 fC± 
an(0) 7 K x \- la lK x \C 2 



(D.53) 



+ 1 
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Lemma |4"31 and the already proven estimates for Vjj+i and vf lead to 



vt 3+1 (z)\ < C 2 r (Cfc(i + 3) + N) 2(i+j+1) i^ 2 +1 Exp(z 



'(i + J + l) 

|021(0)I ff Ci^-2 



a u (0) C 



2e a 7r^2 +C a 1 + 



O2l(0)|\ 7a 1 



+ c 2 \k 2 ) \kJ + 4 c 2 U 2 y U 2 

where the expression in square brackets is bounded by 1 due to our requirement (|D.19|) . 
Next we look at 

l~.L \ 

1 



au(0) J K x l- la /K x \C 2 
~Ka n j+1 



~ L _ ai2(0) ~L , 

an(0) + 



au(0) ai 



iv ; fe=l 



+1 + ag ) (0)t^_ fcii+1 



Again, we proceed by induction on i, the case i = having been handled already. From Lemma [4, 71 
and the induction hypotheses we get 



KVi^l < c x k\. 



' ^i +1 (C 5 (i + j + 1) + |*|) 2 < i+ ' +1 > Exp(z) 



|oi 2 (0)| C 2 , 2e a C a 1 7 a 

+ i — + 



1 + 



G 2 



an(0) C x 1 au(0) an(0) l- la /K x K x \ d. 
Again, the expression in square brackets is bounded by 1 in view of our requirement (|D.2ip . 

Finally, we consider j +1 which satisfies (|4.30p . The right-hand side of the boundary value problem 
satisfies 

min{ij+l} 

RHS KME\ \ < Ca\v^ j+X {z)\+C a ^ lalzf {\utk, j+ l-k\ + \vi-k, j+ l-k\} 

k=l 

< C 3 ^ +1 7r ^ T (C s (i + j) + \z\f^) +1 Exp(z) 



C a 



C 9±k~ 2 + 

"<".i 4 K 3 K 3 l- la /{K 3 K 3 ) 



Lemma 14.51 together with the induction hypotheses, therefore gives us for the solution ufj +1 of 



\u± +x (z)\ < C 3 KlKi +1 ' 



(i + j + 1)1 



(Cs(i + j + l) + \z\ 



\2(i+j+l) 



— + 



Lan(O) i C 3 K 3 K 3 l- la /(K 3 K 3 ) 

+ 9i (EX f J_Y +1 + 9i (EiX (I± 

C 3 \K 3 ) \K 3 ) C 3 \K 3 ) \K 3 



Exp(z) 



By our requirement (|D.22p . the expression in square brackets is bounded by 1. 
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D.4 Proof of Theorem I4A31 



We first prove two auxiliary lemmas. 
Lemma D.l. For every 7 > the functions 

h(k) ■■= ^(i-kf- k (j-ky- k , 
f2(k) ■.= k k 7 k (i-ky- k (j-ky~ k , 

are convex on (0, min{i, j}). 

Proof. It is easy to check that 4bz In f\(k) and ^ln/ 2 (A;) > 0. □ 
Lemma D.2. Let a < b < c < d. Let f be non-negative and convex on [a,d]. Then 

II/IIl°°(6, c ) = max{min{/(a), /(&)}, min{/(c), /(d)}}. 

Proof. We restrict our attention to the case a<b<c<d - the general case can be proved using 
similar arguments. By convexity, we have ||/||i,°°(& jC ) = max{/(6), /(c)}. We claim that 

max{min{/(a), /(&)}, min{/(c), f(d)}} = max{/(6), /(c)}. (D.54) 

Suppose max{min{/(a), f(b)}, min{/(c), f{d)}} < max{/(6), /(c)}. Then 

min{/(a), f(b)} < max{/(6), /(c)} and min{/(c), f(d)} < max{/(6), /(c)}. (D.55) 

If /(a) < f(b), then we write b = Xa + (1 — A)c for some A G (0, 1), and use convexity to get 

/(&) = /(Aa + (1 - A)c) < A/(a) + (1 - A)/(c) < A/(6) + (1 - A)/(c), 

from which we conclude f(b) < /(c). Thus 

/(a) < f(b) < /(c). 

The second condition in (|D.55P then produces min{/(c), f(d)} < /(c), from which we get f{d) < 
/(c). On the other hand, convexity implies upon writing c = Aa + (1 — A)d for some A G (0, 1), 

/(c) = /(Aa + (1 - Ad)) < A/(a) + (1 - A)/(d) < A/(c) + (1 - A)/(c) = /(c), 

which is a contradiction. We conclude that /(a) > /(6). 

Next, we investigate the possibility /(d) < /(c). We proceed analogously. Convexity (for the points 
b, c, d) implies 

f(d) < /(c) < /(&). 

The first condition in (|D.55|) then produces min{/(a), /(&)} < /(&) from which we get /(a) < f{b). 
Using again convexity for the points a, b, c, yields f(b) < /(c), which is the desired contradiction. 
We conclude that /(c) < /(d). 

Combining the above, have /(a) > f{b) and /(c) < /(d) and therefore (ID. 541) . □ 



54 



Proof of Theorem I4.13t We have to estimate the terms on the right-hand side of (|4.40|) . A 

direct check shows that the estimate is valid for the special case M 2 = 0. Hence, we will assume 
M 2 > 1. 

Throughout the proof, we will use that x is real and strictly positive. 
First, we estimate the double sum 



Mi Mi 



j=0 j=M 2 -l 



e ^ 



Using the bounds of Theorem 14.91 and the Cauchy integral theorem for derivatives (with contour 
dBi(x)), we get with the aid of Lemma FOl 



\u'i j {x)\<cpj i+ h\j + iy + \ 

for a suitable constant 7. Therefore, if jJrj{Mi + 1) < 1/2, then 

Mi M 2 /oX j / N M 2 -l 



Mi M 2 . , 

E E m'( £ ) Kii<^( £ 



i=0 j=M 2 -l 

Hence, the double sum in (|4.40p can be estimated in the desired fashion by requiring 7 > 7/2. 

We now turn to the triple sum in (|4.40j) . We start by writing the conditions (|4.39p on the indices 
i, j, and k appearing in the triple sum (14.40p in a more compact form, by setting 

X := {(i,j,k): (i > Mi + 1 V j > M 2 - 1) A max{z - M u j - M 2 } <k< min{i,j}}. (D.56) 

Thus, the triple sum can be written as 

S= £ n\e/»yx k A k ( f^~ k ) . (D.57) 
(i,j,k)ex v u 7 

Using (041 . (g35|) and ||A fc || € i < C A j k , where C A ■= 2C a , we obtain with K > max{i^3, K4}, 
K > max{i^3, K^}, C > max{Cg, Cg} for x > 0, with the aid of Lemma 14.81 



t fr M 



/- \2(i-k+j-k) 

^ ■ u [C(i-k + j-k) + \x\) 

< cc A v m^k^k'-^ pn+jr^i 

< CC A fi i (e/fi) j K i K j ^^Y- k +^- k (i-ky- k (j-ky- k e~ 3 ^ 4 



here, we selected 7 > suitable in dependence on C and a. It is convenient to abbreviate 

k) = ^- k {e^y- k k\- k ki- k {xe la ) k {i - ky- k u - ky- k , 
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with K\ := K'y and K 2 := Kj. Hence, we wish to estimate 

g -3ax/4 f(i,j,k). (D.58) 

(i,j,k)ex 

Next, in order unify the presentation, we consider the cases j = M 2 — 1 and j = M 2 separately. 
That is, we write 

X C XUli, 

Zi := {(i, j, k): j e {M 2 — 1, M 2 } A max{i — Mi,Q} < k < mm{i, j}}, 
X := {{ij, k) : (i > M x + 1 V j > M 2 + 1) A max{i - Mi, j - M 2 } < k < min{i, j}}, 

and estimate the sums over X\ and X separately. 

The structure of the remainder of the proof is as follows: 

• In Step 1, we estimate e~- x ^ Y^(i,j,k)eXi 

• In Step 2, we estimate e"^' /4 J2(ij,k)ei k )> 

• Finally, in Step 3, we consider the case xe > c and show that then \L £tfJ TJg L (x)\ < Ce~^ x . 
Step 1: We estimate 

Mi M 2 min{i,j} oo M 2 min{i,j} 

E /(j lijfc ) = £ £ £ E E E + 

(ij',fc)GXi i=0j=Af 2 -l fc=0 i=Mi + l j=M 2 ~l k=i-M 1 

using the convexity properties of the function /. Specifically, in order to estimate Si, we first 
consider the case M x < M 2 — 1. 

Step la: Assume Mi < M 2 — 1. Then by convexity of the function k 1— > f(i,j, k) (cf. Lemma |D.1|) . 

Mi M 2 i Mi M 2 

5l = E E E/^'^ fc )^E E (i + l)max{/(i,j,0),/(i > i,i)} 

i=0 j=M 2 -l k=0 i=0 j=M 2 -l 

Mi A/ 2 

< (Mi + 1)E E [^\e/^K\K{ef + {e/^' i Ki-\xe la r{j-iy~ l 

i=0 j=M 2 -l 

< C(Mi + l)(e/nK 2 (M 2 - l))^- 1 + C(Mi + l) 2 max{(£ 2£ / M (M 2 - l))^" 1 , (xe 7a ) Ml }, 

where we employed the assumption that \iM\ and e/fj,(M 2 — 1) are sufficiently small and the 
convexity of the second term (as a function of i). From Lemma lD.21 we can bound 

for suitable 7, so that we obtain together with Mi < M 2 — 1 and the trivial bound e < fi, 



e -ax/i Si < c 



f \M2-l 

(£(M 2 -1) 7 J +( M M l7 ) M i 



(D.59) 
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for suitable constants C, 7 > 0. 

Step lb: Analogous reasoning covers the case Mi = M 2 . More precisely, we write in this case 

Mi M 2 min{i,j} Mi-1 M 2 i M 2 j 

E E E E E E + E E/(M,j,fc)=:*,i+^. 

»=0j=Afa-l fc=0 i=0 j=M 2 -lk=0 j=M 2 -ik=0 

For 5i 5 i, the reasoning of the above case ("Step la") is applicable, since Mi — 1 < M2 — 1, and 
yields 



e -ax/A Si i < c 



/I 



M2-1 



(M 2 - 1)7 + (e(Mi - l) 7 ) AJl 



Mi-1 



< C ( -(M 2 - 1)7 



m 2 -i 



where in the last step, we used the trivial bound e < e/jj, in view of fi < 1 and the fact that 
Mi = M 2 . For 5i,2, we use convexity of / in the third argument to arrive at 

M 2 

Si, 2 < (M 2 + l) E max{/(M 1 ,i,0),/(Mi,j,i)}. 

j=M 2 -l 



Estimating max{f(M 1 ,j,0)J(M u j,j)} < f(M 1 ,j,0) + f(M 1 ,j,j), we get 

f(M u j,0) + f(M 1 ,j,j) < ^(e/ t iyK^^M^f+ij, M ^K^- j (xe 7a y(M 1 - 
Using the fact that \xM\ and ej [iM 2 are sufficiently small and that Mi = M2, we get 

M 2 

max{/(Mi,j,0),/(M 1 ,j,j)} < 

j=M 2 -l 

M2-I 



( fJ ,K 1 M 1 ) M ^ ( -K 2 (M 2 - 1) ) + (& 7 a) JUl + /x#i(s£7a 



\M2-1 



Upon writing 



,M 2 -1 



^-7a 



M 2 -l 



/U M2 < ( x- 7a 



M 2 -l 



we can use ()D.59|) to argue as in Step la. 



Step 1c: We now consider S\ for the case Mi > M2 + 1. We write 

M 2 -l M 2 i Mi M 2 j 

*=E E E/fc**)+E E E/fc**)- 

i=0 j=A/ 2 -lfc=0 i=M 2 j=M 2 -l k=0 

Checking the arguments for the case Mi < M2 — 1 of Step la, we see that we can bound 



M 2 -l M 2 i 

e^E E E/fc**)^ G 

i=0 3= M 2 -l fc=0 

which has the desired form in view of [i < 1. 



m 2 -i 

-(M 2 - 1)7 ) + (e(M 2 - 1) 7 ) M2 - 1 
A 1 



57 



It therefore remains to estimate Yli=M 2 X/j=M 2 -i Yli=o /(*>i> which we do again by exploiting 
convexity of the function k \-t f{i,j, k). We get 

Mi M 2 j Mi M 2 

E E E/(*>**0 ^ E E (j + l)max{/(i,i,0),/(i,i,j)} 

i=M 2 i=M 2 -l fc=0 i=M 2 j=M 2 -l 

Mi M 2 

< (Af a + 1)5] ^ ji' /i)' k:{i'j J ■ /,' 'A'i ./;' ' 

i=M 2 j=M 2 — 1 

< C(M 2 + 1) [{e/n{M 2 - W*- 1 + M (M! + l)(xe 7a ) M2 - 1 + (& 7 a) M2 ] 

< C(M 2 + 1) [(e//i(M 2 - lh)^- 1 + (3e7a) Ma_1 ] , 

where we exploited again the assumption that n{M\ + 1) and ej ' \i{M 2 + 1) are sufficiently small so 
that sums can be estimated by convergent geometric series. The contribution M 2 (xE"f a ) M2 ~ 1 can 
now be estimated as before using (|D.59|) . 

Step Id: We now turn to S 2 and start with assuming M\ + 1 > M 2 . Then S 2 takes the form 

oo M 2 j Mi+M 2 M 2 j 

s *= E E E /(<>**) = E E E 

i=Mi+l j=M 2 -l fc=i-Mi 

Convexity of i— > f(i,j, k) allows us to infer 



i=Mi+l j=M 2 -l k=i-M\ 



max /(t,j,fc) < ^^(x^)'"^^) 

K = !-Ml,...J 

+ /i i -^- i (xe7a) i (*-j) < - i . 



i-i+Mi i^j-i+Mi M M 1( 



^^M; 1 (j -i + Mi> 



j-i+Mi 



Since < j — i + M\ < M 2 , we can exploit that e/fi(M 2 + 1) is sufficiently small and since 
< i — j < M\ , we can use that \iM\ is sufficiently small to conclude 



S 2 < C(M 2 + 1) 



Mj+M 2 



Y {liK 1 M 1 ) M ^xe la ) i - M ^ + Y, ( Je 7«) 



i=Mi + l j=M 2 -l 

With Lemma 14.81 we therefore get for suitable 7', since i — M\ < M 2 , 



e- x ' 4 S 2 < C(M 2 + 1) 



Mi+M 2 M 2 

{fiMtK^ Y {zM 2l j- h ' h + Y i eM ^') 
i=Mi+l j=M 2 — 1 

< C(M 2 + 1) {^M l k 1 ) Ah eM 2 + {eM 2l ') M ^ 1 



where, in the last step, we exploited again that \iM\ and eM 2 are sufficiently small. Since M 2 < 
Mi + 1, this last estimate has the desired form. 



Step le: Next, we consider the case M\ + 1 < M 2 . We write 
00 M 2 mm{i,j} 

s * = E E E /(m,*) 

i=Mi+l j'=M 2 -l k=i—Mi 

M 2 — 1 M 2 i 00 M 2 j 

= E E E /(*.**)+ E E E f(i,j,k)-s2,i+sw. 

i=Mi j=M 2 -l k=i-M\ i=M 2 j=M 2 -l fc=i-Mi 
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We note that 

M 1 +M 2 M 2 j 

S ^= E E E 

i=M 2 j=M 2 -l k=i-Mi 

We may bound £2,2 using arguments similar to those of Step Id. The convexity of k 1— > f(i,j,k) 
yields 



u max f(i,j,k) < M Ml i^f 1 (^ £ 7a) ^ - JWl (^)^ ^+iMl ^^ ^+J,2l Mr l (J-^ + M 1 

k=t—M\,...,i 



+ (xs<y a ) 1 . 

Since < j — i + Mi < M\ < M 2 , we can estimate by convexity 



+ 



-K 2 (j-i-Mx) 



j-i-Mi 



< ( 1 + (e/fiK 2 M 2 ) M2 ) < C, 



since we assume that e/ ' \xM 2 is sufficiently small. Hence, we get 

Mi+M 2 

S 2 ,2<C(M 2 + 1) ^K 1 M 1 ) M ^xe la ) i - M ' + (xe la )\ 

i=M 2 

With Lemma B~8l we therefore get 

Mi+M 2 

e-^ /4 S 2 , 2 < C(M 2 + 1) Y, {nKxM^^'ii - Mx)y~ Ml + {er/i)\ 

i=M 2 

for suitable 7' > 0. Recalling that M\ < M 2 , we obtain by assuming that eM 2 is sufficiently small, 
with the aid of geometric series arguments, 

e -9*/4s 2i2 < C{M 2 + 1) lifiKxMx)^ + {ei(M 2 + 1)) M2 



For 52,i we note 
, max f{i,j,k)< 

K=!-Ml,...,l 

(e/fiy- l K$-\j - iy- l (xe la y + ^K^M^ie/^y-'+^ixsjay-^^+^U -i + Mx) 

We note i and j in the definition of the sum 62,1 are such that < j — 1 v iW2 ana iwi ^ j—i-\-ivix ^ 
M2. Hence, this setting simplifies under the assumption that [iM\ and e/[iM 2 are sufficiently small: 

max f(i,j,k) < {xe la y + { l iK 1 M 1 ) Nh {xe la y- M ^{ 

k=i—M\,...,i 

< {xe la y + (^K l M l f h {xe la 

With the aid of Lemma 14.81 we estimate with suitable 7' > , 

m 2 -i 

e -™/±S 2 ,x < {Mx + ^Y.ieiiy + ipKxMxf^ieiii-Mx))^ 11 

i=Mi 



j-i+Mi _ 

< j-i < M 2 and Mi < j-i+Mx < 
and e ' 

^- Ml( s/nM 2 K 2 ) m 
,i-Mi 



i=Mi 

< C(Mi + 1) [(eM 27 ') Ml + (/x^iMi) 
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where we employed again suitable geometric series arguments. Using e = (e/ fj,)fj,, we get 

Mi 



e -s£/4£ 21 < C(M 1 + l) 



ft 



Mi 



-M 2 y 



\Mi 



which has the desired form since e/(j,M 2 is assumed to be sufficiently small. 

Step 2: Before preceding, we point out that we make the assumption 

1 

XE^a < -, 



(D.60) 



as the converse case is covered in Step 3 below. We estimate the contribution arising from the sum 
over I and show 



-ax /A 



E f(ho,k)<c 



— (eM 2 j) M2 + fi(fj,Mu) Ml 



in fact, the reasoning below shows a slightly sharper estimate. 



(D.61) 



Noting the appearance of max{z — Mi , j — M 2 }, we have to study the cases i — Mi > j — M 2 and the 
reverse case i — Mi < j — M2 separately. We restrict our attention here to the case i — Mi < j — M2 , 
since the reverse case is easily obtained using the same arguments (effectively, Mi and M 2 and /i 
and e/fj, reverse their roles). In this situation, we consider the following two subcases separately: 

(i > Mi + 1 V j > M 2 + 1) A (i - Mi < j - M 2 ) A (j - M 2 < k < i) A (i < j), (D.62) 
(i > Mi + 1 V j > M 2 + 1) A (i - Mi < j - M 2 ) A (j - M 2 < h < j) A (j < i). (D.63) 

A key ingredient of our proof is the convexity assertion given in Lemma ID. 11 We recall that a 
non-negative convex function attains its maximum at the boundary (i.e., in the univariate case, at 
the endpoints of an interval). 

Step 2a: We consider the case (|D.62|) . which can be further subdivided into the cases i > Mi + 1 
and j > M 2 + 1. 

Step 2al: We consider the case (]D.62p with the further assumption i > Mi + 1. We get 

(i > Mi + 1) A (i < j) A {i<j-M 2 + Mi) A (j - M 2 < k < i). 
That is, we have to estimate the triple sum 

j <i+M 2 i 

T: = £ £ E /(•".**)> 

i>Mi+l j>i k=j-M 2 
j>i+M 2 -Ah 

where the summation index is i for the outermost sum, j for the middle sum, and k for the innermost 
sum. This triple sum is estimated using convexity of the argument. The innermost sum has at 
most i — j + M 2 + 1 terms and k *— > f(i,j, k) is convex by Lemma lD.ll Hence, we obtain 

j<i+M 2 

T < Yj 12 (i-3 + M 2 + l)max{f(i,j,i)J(i,j,j -M 2 )} 

i>Afi+l j>i 



j>i+M 2 -Mi 



j<i+M 2 



< (Mi + 1) E E max{f(i,j,i),f(i,j,j-M 2 )}, 



i>M 1 +l j>i 

j>i+M 2 -Mi 
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where, in the second step we have used the restrictions on the sum on j to bound i — j + M 2 + 1 < 
Mi + 1. Writing out f(i,j,i) and f(i,j,j — M 2 ), we have 

f(i,j,i) = kt i ( £ /^- i (x S r fa ) i (j-iy- i , 

f(ij,j - M 2 ) = ^+ M %e/ri M *Ki-i +M \xe la y- M *K 2 M *(i - j + M 2 y~^M 2 M \ 

which are again convex functions of j by Lemma ID. II Turning now to the sum on j, we see that 
it has at most M 2 + 1 terms. In view of Lemma lD.2( we bound for the relevant j: 

f(i,j,i) <A:=max{f(i,j,i)\ j=i+ M 2 ,f(hj,i)\j=i} = max{/(i, % + M 2 , i), f(i, i, i)}, (D.64a) 
f(i,j,j ~ M 2 ) < B := max{f(i,j,j - M 2 )\j=i+M 2 ,f(i,j,j - M 2 )\j=i+M 2 -M 1 } 

= max{/(i, i + M 2 ,i),f(i, i + M 2 - M u i- Mi)}. (D.64b) 

More explicitly, these are 

A < {xe la ) i mc^{l,K^ 2 {e/^) M2 M 2 M2 }, 

B < (e//i) M2 M 2 M2 ^ 2 M2 m a x{(xe 7 a) i ,(^iMi) Ml (^7a) i - Ml }- 
Writing out the sum and using convexity of the argument, we can estimate 

j<i+M2 i 

S := (M x + 1) J2 E E ma X {f(iJ,i),f(i,j,j-M 2 )} 

t>Afi+l j>i k=j-M 2 
j>i+M 2 -M 1 

00 

< (Mi + l)(min{Mi,M 2 } + l) ^ A + B. 

i=Mi+l 

Using the facts that /uMi and e/ are sufficiently small, we can estimate with geometric series 
arguments, in view of our assumption (jD.60[l : 

S < C(Mi + l)(min{Mi,M 2 } + 1) [(xe 7 a) Ml+1 + {e/ iiM 2 K 2 ) M *xe la 

Hence, using Lemma 14.81 to control x Ml+l in the first term and x in the second term, gives with 
suitable 7', 



e -a$/4 s < c(Mi + i)( min { Ml) m 2 } + i) \{e[M x + l) 7 ') Ml+1 + e(e/fiK 2 M 2 ) M2 
< C 



(e(Mi + l) 7 ") Ml+1 + £( e /^ 7 "M 2 ) M2 



A* 

where, in the second step, we selected 7" suitable and used the fact that \iM\ can be controlled. 
Since e < fj,, we obtain the desired estimate (|D.6ip . 

Step 2a2: The other case is j > M 2 + 1. This leads to 

(j > M 2 + 1) A (i < j) A (i < j - M 2 + Mi) A (j - M 2 < k < i). 
Writing out the sum, we have 

i<j 

i<j-M 2 +hh i 

r := E E E 

j>A/ 2 +l i>j-M 2 k=j-M 2 
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The innermost sum has at most M 2 + 1 terms which, by convexity, can be estimated by 

max{/(«,j, i),f(i, j,j -M 2 )}. 
Writing out f(i,j,i) and f(i,j,j — M 2 ), we have 

f(ij,i) = {Eiii)^ky\xe la )\i-iy-\ 

f(i,j,j - M 2 ) = ^-3+ M 'K{- j+M '(i - j + M 2 y-^ M ^e/fi) Ak K^(xe la y- M ^M^\ 
which are again convex as functions of i. Hence, by Lemma lD.21 we bound for the relevant j: 

f(i,j,i) < A' :=max{f(i,j,i)\ i= j-M 2 ,f(i,j,i)\i=j} 

= m^{(e/fiM 2 K 2 ) M ^xe 7a y- M ^(xe 7a y}, 
f(h3,j ~ M 2) < B' := max{/(«, j, j - M 2 )\i=j-M 2 , f(i,3, j ~ M 2 )\i=j-M 2 +M 1 } 
< (e/fiM 2 K 2 ) Ah (xe-/ a y- Ah max{l, {^hK^}}. 

The middle sum in T' has at most min{Mi,M2} + 1 terms. Hence, we arrive at 

T' < (M 2 + l)(min{M 1 ,M 2 } + l) {e/nM 2 k 2 ) M \xe la y- M * + {xe la y 

j>M 2 + l 

< C(M 2 + l)(min{Mi,M 2 } + 1) [(e/ ^M 2 k 2 f h (xe la ) + {xe la f h+l ] , 

where, in the second step we have employed geometric sum arguments, which are applicable in view 
of (|D,60|) . Reasoning as at the end of Step 2al, we get for suitable 7', 

e -ax/4 T > < C (M 2 + l)(min{Mi, M 2 } + 1) \s{e / iiM 2 K 2 ) h ' h + (e(M 2 + 1) 7 ') M2+1 
< C [e(e/nM 2l ") M2 + (e{M 2 + 1) 7 ") M2+1 ] , 

where, we employed the assumption that fJ,M\ and ej fiM 2 are sufficiently small and, in the second 
step, 7" is selected appopriately. This leads to the desired estimate (|D.6ip . 

Step 2b: The case (|D.63|) is further subdivided into the cases i > Mi + 1 and j > M 2 + 1. 

Step 2b 1: In the fist case, i > Mi + 1, we get 

(i > Mi + 1) A (j < i) A (i < j — M 2 + Mi) A (j - M 2 < k < j). 
We immediately see that we need 

M 2 < Mi, 

for this set of indices to be non-empty. 

Writing out the sum and using convexity of the argument, we get 

t" ■■= E E E /&;.*)• 

i>Mi+l j>i+M 2 -Mi k=j-M 2 
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As before, the innermost sum has M 2 + 1 terms and convexity yields, for the terms of the innermost 
sum, the upper bound max{/(i, j, j), f(i, j, j — M 2 )}. More explicitly, 



f(hj,j) = l< j K\ J U:-„} J (i j) 
f(i,j,j-M 2 ) 



1 1 



{xe la y- M *(e/v) M *K^M^(i-j + M 2 ) 



M 2 jfM 2 ]\/t M 2 1 



\i-j+M 2 



Again, we recognize the functions to be convex in j, so that for the relevant indices j we have 

f(i,3,j) < A" :=max{f(i,j,j)\ j=i ,f(i,j,j)\ j=i+ M 2 -M 1 } 

< max{(xe 7 a) i , (xe la y +M ^ h {^M x - M 2 )K{) M ^ M ^ 
f(i,j,j ~ M 2 ) < B" := max{/(i, j,j - M 2 )\j=i, f(i,j,j - M 2 )\ j=i+M2 _ Ml } 

< m&x{(fj>M 2 Ki) M2 {xe^/ a ) i ~' M2 (e/ [iM 2 K 2 ) M ' 2 , (jiM\Ki) Ml (xej a ) i ~ Ml (e/ fiM 2 K 2 ) M2 }. 

The middle sum of T" has at most Mi — M 2 + 1 terms. Therefore, we get with our assumption 
that \xM\ and ej jiM 2 are sufficiently small and the assumption (|D.60p . 

T" < (Mx - M 2 + 1)(M 2 + 1) [(xe 7 a) Ml+1 + (x£7a) M2+1 (/i(M! - M 2 )K x f h - M2 

+ {ixM 2 k x ) M2 {xe la ) M ^ x - M2 {elyiM 2 k 2 )^ + (/,M 1 ^ 1 ) Ml (^7a)(e//uM 2 ^ 2 ) M2 

Since M\ > M 2 , we can estimate further with the aid of Lemma 14,81 

e -ax/4 T „ < C ( Mi _ M2 + !)( Ml + X ) \(e(Mi + l)7) Ml+1 



+ (e(M 2 + 1) 7 ') M2+ V(^i - M 2 )K X ) M ^ M2 

+ {eM 2 K 1 K 2 ) M2 {e{M 1 -M 2 + i^'fh-^h+i + e ( /U M 1 K 1 ) Ml (e//iM 2 K 2 ; 



M 2 



c 



( £ (M 1 + l) 7 / ) AJl+i + -(e(M 2 + l) 7 / ) 



i\M 2 



where we used again that \xM\ and e / ^M 2 are sufficiently small. This leads to the desired estimate 
(fLT6T|) . 

^ep 565: The last case is j > M 2 + 1. This reads: 

(j > M 2 + 1) A (j < i) A (i<j- M 2 + Mi) A (j - M 2 < k < j). 
We see again that 

M 2 < M u 

is a necessary condition for the set of indices to be non-trivial. 



Writing out the sum we have 

rp/// 



j-M 2 +hh j 

E E E 

j>M 2 +l i>j k=j-M 2 



Again, we use convexity of k 1— > f(i,j, k) and observe the estimates 

f(ij,j) = ^- j k{^{xe la )\t-jr^ 
f(i,j,j - M 2 ) = ^^+ M ^r J+M2 (^7a) J - M2 (e/^) M2 ^ 2 M2 M 2 M2 (i - j + M 2 y-J+ M2 . 
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Since these functions are convex functions of i, we can estimate for the relevant indices j: 



f(i,j,j) < A'" := max{f(i,j,j)\ i=j , f(i,j,j)\ i=j _M 2+Ml } 

< {xe la y max{l, - M 2 )K 1 ) M ^}, 

f{i,j,j - M 2 ) < B'" := max{f{i,j,j - M 2 )\ i=j J{i,j,j - M 2 )\ l=j _ M2+ M 1 } 

< (xs 7a y- M He/vM 2 K 2 ) M * max{ (fiM 2 K\ ) M ' 2 , (fxMi K\ ) Ml } . 

Furthermore, in the triple sum T'" the number of terms in the innermost sum is M 2 + 1 whereas 
the number of terms of the middle sum is bounded by M\ — M 2 + 1. Hence, we can bound with 
geometric series arguments, in view of the assumption (ID.60j) . 



T'" < C{M 2 + l){M 1 -M 2 + l) 



1M2+1 



{xe la ){e/tiM 2 K 2 ) M * max^^i^)^ (//M^)^} 



Since M 2 < Mi, we can estimate \xM 2 K\ < [iM\K\ < 1, where we appealed in the last step to our 
standing assumption that \iM\ is sufficiently small. Hence, max{([iM 2 Ki) M2 , (nMiKi) Ml } < 1. 
Thus, we can simplify 

T'" < C(M 2 + l)(Mi - M 2 + 1) [(xe 7a ) M2+1 + (xe la )(e/fiM 2 K 2 ) M2 

Hence, with the aid of Lemma 14.81 

e -a$/4 T m < C{M2 + 1)(Mi _ M2 + i) \{ e {M 2 + 1) 7 ') M2+1 + e{e/fiM 2 K 2 ) M2 



< C- 



fj, l 



(e(M 2 + 1) 7 ") M2 + (£//xM 27 ") M2 



This leads to the desired estimate (|D.61|) . 

Step 3: We now cover the case when xe is bounded away from zero. Specifically, let c > be fixed 
and consider the case xe > c. Then we have the pointwise estimate 



Mi M 2 

Uto| < ^J> i (£/^(|2.§(x)| + 

i=0 j=0 



. . _,(<?(* + i) + |£ 

i=o j=o 

Mi M 2 



2(i+j) 



:J 0/ i +i 7 -S-ip-3a$/4 



i=0 j=0 



Mi 



Ai 2 



i=0 
Mi 



j=0 



M 2 



< 



e -3aJ/4 J- ^ KlMi y ((£/n)K 7 M 2 ) 
3=0 



i=0 



from which the desired result follows, provided fiK^Mi < 1 and (e/ fj,)KjM 2 < 1. Completely 
analogously, we get bounds for the derivatives of lJg L ■ 
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In view of the form of the differential operator L £41 applied to functions of x (see (|4.10p ). we get in 
view of [x > e, 

\L e ^ L (x)\ < [e 2 + /j 2 ] e- 2 Ce-^' 4 < C^e^e^l 4 < Ce'V^/ 4 . 

From the assumption xe > c, we see that e~ 2 < c~ 2 x 2 , and the factor x can again be absorbed by 
the exponentially decaying g -3 ^/ 4 . □ 



D.5 Proof of Theorem 14.141 



We first study the case 



< xfi < 1/2. 



(D.65) 



The starting point is the expression for L £ ^XJ^ L in (|4.47[) . It consists of a double sum and a triple 
sum. We first consider the double sum. From the bounds on Ujj of Theorem 14.91 and Cauchy's 
integral theorem for derivatives as well as Lemma 14.81 we get for suitable 7' > 0, 

Mi M 2 +2 Mi 

£ £ fl \e/tf\tiZ j _ 2 (x)\ < Ce-^ 2 ^2S(i + M 2 + l) 7 'y +M2 (e/v) M2+1 

i=0 j=M 2 +l i=0 

< Ce-^(K£/n(M 2 + 1)) M2+1 , 
where, in the second step we used Lemma 14.81 again. 

We next turn to the triple sum in (|4.47j) . From (|4.32|) . (|4.33|) . Lemma l4~8l and ||A fc || < Cai\, we 
obtain with K > maxji^i, -f^}) > max{i£i, ^2}) C > max{C^, C%}- 



M 



(i-fc+i)! 



i=Mi+l j'=0 fe=i-Mi 
00 M2 i 

' rr > E E E S(e/ri j ^~ k K>x k f + i- k (i + j~k) 

i=Mx+l j=0 k=i-Mi 



2(i-k+j) 



-aRe(i) 



i+j— kg— a3x/A 



The argument is a convex function of k. Hence, we can bound for suitable 7 > 0, 



T TT M 



00 M 2 



< Ce~^ x/ \M l + l) Y^{s/fj,y(K^yU j x H + x i - Ml K Ml {j + M 1 y +Ml 

i=Mi+lj=0 



In view of our assumption xfj, < 1/2, the outer summation on i leads to a convergent geometric 
series. For the inner summation, we use (j + Mi)i +Ml < p M^ x e^ Ml ^ the assumptions that \iM\ 
and e/ ' [1M2 are sufficiently small and convergent geometric series arguments to get, with appropriate 

7>0, 



T TT M 



< Ce-^ x/i {M l + 1) [(fix) Ml+1 + (/zMn) 7 
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which can be estimated in the desired fashion with the aid of Lemma 14,8 



T TT M 



< Ce~^ /2 (M l + l) e -^'/4 [{^x) Ml+1 + (^Mi7) Ml Fj 

< Ce-^l 2 [(/i 7 '(Mi + l)) Ml+1 + (/t(Mi + l)7) Ml+1 ] , 
which has the appropriate form. 

We now consider the converse case xfi > 1/2. Here, we need to exploit the fact that the functions 
U|f L are exponentially decaying. We observe that the same reasoning as in Step 3 of the proof of 
Theorem 14.131 yields 



< Ce" 3 ^ /4 , 



and, by Cauchy's integral theorem for derivatives, estimates for the derivatives. The arguments of 
Step 3 of the proof of Theorem 14.131 therefore, yield 
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